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Figure  1.1  Three  views  of  the  neutral  curve  for  an  example  parallel  shear  How:  (a) 
temporal  stability  approach,  (b)  spatial  stability  approach  and  (c)  spatial  bifurcati 
approach. 


Figure  1.2  Effect  of  the  wall  elastic  modulus  on  the  neutral  curve  for  the  Blasius 
boundary  layer  adjacent  to  a  compliant  wall  (after  Carpenter  &  Garrad  (1985,  Figure 
11]). 


Figure  2.1  Neutral  curve  in  the  (c,il)  plane  for  the  (parallel)  Blasius  boundary  layer. 

Figure  2.2  Spectrum  of  the  modified  (real)  Orr-Sommerfeld  equation  for  fixed  c  € 
(ci,C2)  as  R  intersects  the  neutral  curve. 

Figure  2.3  Possible  movement  of  the  spatial  Floquet  multipliers  exp(2sr7/a)  as  a 
function  of  e. 


Figure  2.4  Global  loop  structure  for  fixed  c  of  spatially  periodic  states  bifurcating 
from  the  Blasius  boundary  layer:  (a)  supercritical  loop  and  (b)  subcritical  loop. 

Figure  2.5  Global  period-doubling  loops  with  a  finite  cascade  in  the  map  (2.15):  (a) 
7  <  1  resulting  in  an  absence  of  period-doubling  and  (b)  7  >  1  (in  particular  7  =  1.30) 
resulting  in  three  period-doubing  bifurcations. 

Figure  2.6  Finite  and  infinite  period-doubling  cascade  in  the  map  (2.16)  where  the 
primary  loop  is  subcritical  with  m  =  /?  =  —y/y  —  .2  and  (a)  7  =  .21  and  (b)  7  =  .24. 

Figure  3.1  Neutral  curves  in  the  (c,  R)  plane  for  the  modified  (real)  3D  Orr-Sommerfeld 
equation  (3.7)  for  /?  >  0. 


Figure  4.1  Neutral  curve  of  the  Orr-Sommerfeld  equation  for  ft  =  0  and  0^0 


illustrating  the  codimension  2  intersection  point. 
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Figure  5.1  Effect  of  reducing  the  wall  clastic  modulus  ( E )  on  the  neutral  curve  for 
the  Blasius  boundary  layer  (after  Carpenter  &  Garrad  [1985,  Figure  11]). 

Figure  5.2  Neutral  curve  at  the  critical  value  of  the  wall  elastic  modulus  E  =  E0  in 
the  (c,  iZ)  plane. 

Figure  5.3  Schematic  bifurcation  diagrams  for  the  normal  form  in  equation  (5.8) 
showing  how  secondary  bifurcations  to  2-tori  and  3-tori  arise:  (a)  infinite  branch  of  T2 
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Figure  5.4  Coalescence  of  the  non- resonant  Hopf-Hopf  interaction  by  the  addition  of 
a  third  parameter  producing  a  1  : 1  non-semisimple  Hopf. 


1.  Introduction 


Many  of  the  classic  fluid  flows  of  great  practical  importance:  boundary  layers, 
channel  flows,  jets  and  wakes  arc  open  systems.  That  is,  they  arc  unbounded  in  some 
(or  all)  spatial  directions  and  therefore  do  not  appear  to  be  natural  candidates  for 
application  of  finite  dimensional  dynamical  systems  theory.  On  the  other  hand  the 
linear  theory  of  such  systems  is  well  understood  particularity  when  the  initial  instability 
is  connective  rather  than  absolute.  The  useful  feature  of  convectively  unstable  flows  is 
that  changes  ill  the  flowfield  take  place  in  a  “boosted”  frame  of  reference  x  t->  x  —  ct 
c  €  H.  The  idea  of  convective  instability  is  intimately  linked  with  the  now  classic  theory 
of  spatial  stability  of  shear  flows  (Caster  [19G3]).  We  introduce  a  natural  generalization 
(to  the  nonlinear  regime)  of  the  theory  of  spatial  stability  (spatial  bifurcation  theory) 
and  our  claim  is  that  the  transition  to  turbulence  in  open  systems  with  equilibrium 
slates  initially  unstable  through  a  convective  instability  can  be  analyzed  using  dynamical 
systems  theory  in  a  spatial  setting. 

In  the  linear  theory  of  spatial  stability  the  frequency  u  of  a  disturbance  is  treated 
as  real  and  the  (in  genend  complex)  wavenumber  is  the  eigenvalue  (with  Im(o)  <  0 
corresponding  to  a  spatially  unstable  wave  and  Im(ar)  >  0  corresponding  to  a  spatially 
stable  wave).  In  the  neutral  case  the  temporal  and  spatial  linear  theories  coincide  and 
a  neutral  curve  can  be  plotted  in  three  ways  as  shown  in  Figure  1.1.  Figure  1.1(c)  is  the 
relevant  neutral  curve  for  spatial  bifurcation  theory  however.  In  particular,  in  spatial 
bifurcation,  theory  the  vtavespeed  c  is  treated  as  a  given  real  parameter.  A  sketch  of  the 
theory  is  as  follows.  Suppose  ( U(y),[ ))  is  a  (parallel)  2D  equilibrium  stale  and  write 
the  Navior-Slokes  equations  (in  this  case  the  stream  function  and  vorlicily  variables) 
perturbed  about  the  equilibrium  state  as  a  spatial  evolution  conation  in  the  boosted 
frame  x  >■>  a:  —  ct, 
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or  succinctly, 

r\ 

|U‘  =  L(c,/i)  •*  +  •••  (1.2) 

when*  t>  -  £  is  the  vorticity  and  w  =  (see  Section  2.1  for  dot  tills).  Tin;  idea  l 

is  to  pick  any  c  £  R  (interesting  values  arc  those  that  intersect  the  neutral  curve), 
increase  R  and  determine  all  bounded  (in  the  streatmvise  direction)  solutions.  In  this 
setting  it  is  straightforward  to  use  the  symmetry  of  the  equations  and  to  apply  dy¬ 
namical  systems  theory  to  show  the  existence  of  a  spectacular  zoo  of  spatial  struc¬ 
tures  including  spatially  quasi-periodic  states  with  2,3  and  possibly  4  independent 
wavenumbers!  The  bifurcation  sequence  begins  at  the  neutral  point  R.  =  R„  where 
L(«,  /?„)<]•(»/)  =  iaA'{y)  ««  <5  R;  that  is,  the  linearization  of  (1.1)  luus  purely  imaginary 
eigenvalues,  a  two-dimensional  center  subspace  and  a  point  of  spatial  llopf  bifurcation. 
With  no  restriction  (except  for  boundedness)  placed  on  the  streatmvise  spatial  structure 
the  spatially  periodic  state  will  inevitably  undergo  wavelength  doubling  (with  cascades) 
and/or  secondary  bifurcation  to  spatially  quasi-periodic  states.  In  fact  a  central  obser¬ 
vation  of  our  work  is  that  spatially  quasi-periodic  states  are  prevalent  in  shear  Hows 
(genetically  occur  in  the  one  parameter  family  of  2D  states  along  the  upper  branch  of 
the  2D  neutral  curve  in  Figure  1.1(c)).  Although  the  theory  presented  in  the  sequel 
is  generally  applicable  to  any  2D  parallel  equilibrium  state  with  a  neutral  curve  as 
in  Figure  1.1  we  suppose  throughout  that  the  basic  equilibrium  slate  is  the  (parallel) 
Ulasius  boundary  layer. 

Of  fundamental  importance  in  transitional  shear  Hows  is  the  origin  and  subsequent 
bifurcation  of  3D  states.  However  in  Section  2  we  begin  with  spatial  bifurcations  in 
2D  and  show  that  even  in  2D  new  and  interesting  spatial  structure  arises.  The  2D 
Navier-Slokes  equations  can  be  written  as  a  spatial  evolution  equation  in  a  number  of 
ways  and  two  forms  arc  introduced  in  Sections  2.1  and  2.3  using  the  stream  function 
&  vorticity  variables  and  the  primitive  variables  (which  leads  to  an  interesting  non- 
standard  evolution  equation)  respectively.  The  basic  2D  spatial  bifurcation  problem  is 
introduced  in  Section  2.1  and  in  Section  2.2  a  spatial  secondary  “instability”  theory 
is  introduced  that  complements  the  temporal  secondary  instability  theory  of  Orszag  &  ' 
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Patera  [1083]  and  Herbert.  [1983, 198-1].  It  is  shown  how  the  known  slructmc  of  wave¬ 
length  doubling  will  potentially  lead  to  cascades  of  wavelength  doubling  (wavelength 
“bubbling'’)  and  the  secondary  bifurcation  to  2D  spatially  quasi-periodic  states  is  ex¬ 
pected  (a  demonstation  of  secondary  bifurcation  to  2D  spatially  quasi-periodic  states 
is  carried  out  in  Section  5). 

Spatial  bifurcations  with  the  addition  of  spamvisc  structure  (three  dimensionality) 
arc  considered  in  Section  3.  The  3D  Navicr-Stokcs  equations  arc  written  as  an  evolution 
equation  in  the  primitive  variables;  that  is  with  'I‘  =  (u,  t>x,  to,tor,p,j>r)y  the  3D  Navicr- 
Stokos  equations  can  be  written  as  =  L(c,/I)'P  +  N($,u;/l)  and  u  is  obtained 
from  the  streamwise  momentum  equations  (see  Section  3.5  for  details).  We  have  not 
attempted  to  construct  other  (spatial)  evolution  equations  for  the  31)  Navicr-Stokes 
equations  but  spatial  evolution  equations  for  the  vorticity  &  velocity  formulation  or 
a  vector  stream  function  formulation  should  also  be  useful.  Any  bounded  spamvisc 
structure  (satisfying  the  equations)  is  admissablc  but  with  the  simple  assumption  of 
spamvisc  periodicity  already  the  number  of  interesting  bifurcations  or  the  stjeamwisc 
structure  is  considerable.  The  assumption  of  spamvisc  periodicity  leads  to  an  0(2)- 
equivariance  of  the  evolution  equation  which  is  central  to  the  analysis  of  bifurcating 
3D  states.  In  Sections  3.1  to  3.3  0(2)-equivariant  (spatial)  Ilopr  bifurcation  theory 
is  used  to  analyze  the  primary  and  secondary  spatial  bifurcation  of  31)  states  that 
arc  periodic  in  both  the  spamvisc  and  streamwise  directions.  Section  3.3  contains  a 
generalization  of  the  spatial  secondary  “instability”  theory  of  Section  2.2  to  3D.  Our 
most  useful  observation  with  regard  to  applications  is  that  all  along  the  upper  branch 
of  the  2D  neutral  curve  there  exists  an  interaction  between  a  2D  state  with  stieaimvise 
wavenumber  cvj  and  a  3D  state  with  streamwise  wavenumber  02  but  with  both  waves 
travelling  at  the  same  phase  speed.  Fr  om  a  theoretical  point  of  view  the  interaction  is 
a  codimension  2  point,  of  an  0(2)-cquivariant  vectorfield  with  a  G-dimensional  center 
subspace!  I11  Section  4  a  formal  application  of  centre-manifold  theory  and  normal  form 
theory  is  used  to  show  that  all  along  the  upper  branch  of  the  2D  neutral  curve  there 
exists  secondary  bifurcation  to  3D  states  that  are  quasi-periodic  in  the  streamwise 
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direct  ion  (and  periodic  in  t  he  spamvise  direction).  The  theory  is  formal  simply  because 
the  Blasius  boundary  layer  is  not  an  exact  solution  of  the  Navier-Slokes  equations 
and  the  additional  neglect  of  non-parallel  terms.  For  strictly  parallel  Hows  with  a 
similar  neutral  curve  (such  as  plane  Poiseuillc  flow)  the  theory  can  be,  carried  through 
rigorously  (Bridges  [1991c]),  although  particular  care  is  always  necessary  when  dealing 
with  the  bifurcation  of  tori. 

The  theory  for  the  quasi-pcriodic  interaction  of  a  2D  wave  with  2  oblique  (3D) 
waves  is  of  gieat  practical  interest  because  it  is  a  mathematically  consistent  theory  for 
the  appearance  of  quasi-pcriodic  waves  in  the  Blasitts  boundary  layer.  Experiments  of 
Kachanov  &  Levchenko  [1984]  have  shown  that  a  quasi-pcriodic  interaction  between  a 
2D  fundamental  Tollmien-Schlichting  wave  with  a  pair  of  oblique  waves  is  observed  as  a 
lobtist.  pat  t  of  th"  transition  process.  The  normal  form  for  the  quasi-pcriodic  interaction 
is  worked  out  in  Section  4.  Some  straightforward  (although  lengthy)  calculations  tire 
necessary  to  determine  the  coefficients  in  the  normal  form  and  this  work  is  in  progress 
(Bridges  [199 lb]). 

In  Section  4.2  the  interesting  idea  of  spanwisc  resonances  is  considered  briefly.  In 
othci  words  two  pairs  of  oblique  waves,  one  vnl:  rpanwise  wavenumber  ft  and  other 
with  span  wise  wavenumber  nft  #  =  2,3,...  interact.  This  is  a  codimension  2  interaction 
(plot  the  ft  and  nft  neutral  curves;  the  point  of  intersection  between  the  two  curves  is  the 
interaction  point).  Such  a  codimension  2  point  occurs  for  each  value  of  ft  ( ft  sufficiently 
small)  but  the  strcamwisc  wavenumbers  of  the  two  waves  will  be  irrationally  related. 
Although  the  spamvise  resonant  interactions  occur  at  Reynold’s  numbers  considerably 
higher  than  the  2D-3D  interaction  of  Section  4.1,  they  are  nevertheless  of  great  interest. 
From  a  theoretical  point  of  view  the  interaction  corresponds  to  an  8-dimensional  ccntre- 
subspaee  and  the  normal  form  indicates  the  potential  for  bifurcation  to  high  dimensional 
tori.  From  a  practical  point  of  view  the  spamvise  resonances  introduce  new  spatial 
structure  Unit  may  lie  important  for  the  transitional  boundary  layer  at  higher  lleynolds 
number. 

Finally  in  Section  5  the  strictly  two-dimensional  problem  is  reconsidered  and  the 
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“codimcnsion-2  strategy”  is  used  to  show  secondary  bifurcation  to  2D  spatially  quasi- 
pcriodic  slates.  Tlic  compliant  wall  problem  (Carpenter  (1990),  Carpenter  &  Morris 
[1990))  pi  ovules  an  interesting  setting  for  the  analysis  because  it  already  contains  nu¬ 
merous  new  parameters.  Research  of  Carpenter  &  Garrad  [1985]  has  shown  that  the 
upper  and  lower  branches  of  the  neutral  curve  coalesce  and  detach  (see  Figure  1.2)  when 
the  elastic  modulus  of  the  compliant,  wall  (adjacent  to  the  Blasius  boundary  layer)  is 
reduced.  In  Section  5  the  critical  point  E  =  Ea  is  analyzed  and  it  is  shown  that  at 
the  interaction  point  of  the  upper/lowcr  branches  the  linear  Navicr-Stokes  equations 
have  spatially  quasi-periodic  solutions  with  two  independent  wavenumbers.  Applica¬ 
tion  of  dynamical  systems  theory  (in  particular  Section  7.5  of  Guckenheimer  &  Holmes 
[19S3|)  shows  that  the  unfolding  of  the  above  singularity  results  in  a  secondary  bifur¬ 
cation  of  (spatially)  quasi-periodic  slates  (in  the  nonlinear  equations).  The  singularity 
in  question  is  not  particularly  important  to  the  main  function  of  the  compliant  wall 
(stabilization  and  drag  reduction)  but  it  nevertheless  demonstrates  that  nr.condu.ry  bi¬ 
furcation  to  spatially  quasi-periodic  states  is  to  be  anticipated  even  in  two-dimensional 
hound  ary  layers. 

In  spite  of  the  fact  that  our  methods  are  local  (centre-manifold  theory,  normal 
form  theory  and  local  cqui variant  dynamical  systems  theory)  the  existence  of  quite 
complex  spatial  structures  in  shear  Hows  is  demonstrable.  On  the  other  hand  il  is  clear 
that  the  bifurcations  to  the  various  spatial  tori  and  sequences  of  wavelength  doubling 
will  inevitably  lead  (in  some  regions  of  parameter  space)  to  chaotic  spatial  structure. 
Will  this  lie  related  to  turbulence?  Suppose  that  a  sequence  of  bifurcations  takes 
place  leading  to  non-trivial  spamvise  variation  and  chaotic  streamwise  structure.  The 
llowlield  will  indeed  be  three-dimensional  and  although  the  governing  equations  aie 
“steady”  the  streamwise  coordinate  is  in  fact  x  —  ct.  Therefore  a  probe  at  lixed  x 
will  legist cr  a  chaotic  How  in  time  even  though  time  does  not  appear  independently. 
Assuming  there  is  suilicieut  three-dimensionality  for  true  turbulence  (and  an  absolute 
instability  does  not  occur)  it.  is  cntiicly  possible  that  the  transition  process  takes  place 
in  the  convective  frame  x  —  ct.  We  call  this  structure  convective  chaos  or  convective 


turbulence  if  indeed  it.  is  turbulence.  However,  in  studying  the  sequence  of  bifurcations 
in  Hie  convective  frame  it  is  important  to  check  for  secondary,  tertiary,  etc.  absolute 
instabilities  which  will  force  time  into  the  problem  independent  of  *  —  et. 

With  the  emphasis  throughout  on  spatial  bifurcations  and  spatial  invariant  man¬ 
ifolds  the  stability  assignments  obtained  from  the  normal  forms  (written  as  evolution 
equations  in  x)  will  not  be  applicable.  A  complete  stability  analysis  of  the  spatial 
struct  ures  is  a.  very  interesting  problem  and  will  require  the  rcintroduclion  of  time  and 
consideration  of  sideband  instability  (and  its  generalizations).  For  spatial  states  that 
are  periodic  in  both  the  stroamwiso  and  spanwisc  direction  the  spatio-temporal  sta¬ 
bility  can  be  studied  using  the  theory  of  sideband  instability  but  the  generalization  of 
this  concept  to  study  the  stability  of  the  spatially  qnasi-pcriodic  states  is  by  no  means 
clear  and  will  be  an  interesting  area  for  further  research. 
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2.  Spatial  bifurcations  in  two  dimensions 


The  2D  Navier-Stokc®  equations  arc  considered  in  a  steady-frame  (a:  i — »  a;  —  ct) 
with  c  given  and  we  look  for  all  possible  states  moving  at  speed  c  but  with  various  spa¬ 
tial  structure  that  bifurcate  from  the  equilibrium  state.  For  definiteness  suppose  that 
the  basic  equilibrium  state  is  the  Blnsius  boundary  layer  (U,V)  with  IJ(x,y)  —  /'(»/), 
V(x,y)  =  T(a-/i)— (»//'(»/)  -  /),  »/  =  y(R/x)1/2  and  /(»;)  satisfying  the  Blnsius  equa¬ 
tion  2/"'  //"  =  0  on  »;  G  |0,oo).  The  Blnsius  solution  is  a  troublesome  equilibrium 

state  in  that  it  satisfies  the  Navier-Stokes  equations  only  asymptotically  and  a  proper 
existence  theory  and  stability  theory  requires  careful  use  of  asymptotic  and  multiple 
deck  theoiy  (Smith  [1979]).  We  assume  here  that  a:  =  O(Il)  and  invoke  the  parallel  flow 
approximation.  Then  U  —  C?(l),  V  =  0(R~l)  and  the  two-dimensional  Navier-Stokes 
equations  perturbed  about  the  (parallel)  Blnsius  boundary  layer  solution  arc 


£  +  £=  0 

ox  dig 

...  . du  ..  dy  1  .  du.  du 

(U-c)^+V,,+  ^-1Au  +  «-s-  +  „- 

. dv  dp  l  dv  dv 

(u  ”  c)oi  +  ojl _  7i A” +  "SJ  +  "aj  =  °' 


=  0 


(2.1) 


In  the  case  where  the  equilibrium  state  is  strictly  parallel  the  set  (2.1)  is  exact.  In  the 
set  of  equations  (2.1)  time  has  been  eliminated  by  the  shift  x  i->  x  -  ct;  in  particular  we 
are  looking  for  steady  bifurcations.  Time  can  be  reintroduced  for  a  stability  analysis 
or  if  there  is  a  bifurcation  to  a  noil-trivial  temporal  state.  The  crucial  difference  in 
(2.1)  from  the  usual  theory  is  that  c  is  treated  as  a  given  real  number.  The  idea  is  to 
treat  the  set  (2.1)  as  an  evolution  equation  in  a;  that  is,  a  pde  with  y  as  the  “spatial” 
variable  and  x  as  the  “time-like”  variable.  With  parameters  (c, R.)  we  look  for  all  the 
usual  bifurcations  in  evolution  equations:  Ilopf- bifurcation,  period-doubling  and  torus 
bifurcation,  etc,  except  that  in  the  spatial  setting  these  bifurcations  will  correspond  to 
a  spalial-llopt  bifurcation,  wavelength  doubling  and  spatially  quasipcriodic  .states. 

The  basic  bifurcation  from  the  Blnsius  boundary  layer  is  determined  by  studying 
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the  spectrum  of  the  linearization  of  (2.1), 


£“  +  =  o 

Ox  di/ 

_ ‘a.-0 


(2.2) 


Now  let  (u,n,p)  =  cA*(ti, «,/>)  then  with  pyy  +  A2p  =  —2\UyV  the  set  (2.2)  can  bo 
reduced  to  a  modified  (rc.il)  form  of  the  familiar  Orr-Soinmerfcld  equation 

L  •  f>  =  (JL  +  A2)  v  +  XHUyyV  -  RX(U  -  c)  (JL  -1-  A2)  S  =  0  (2.3) 

We  have  purposefully  used  cXr  rather  than  the  usual  c,ttX  (which  we’ll  revert  back  to 
shortly)  (o  make  an  analogy  with  the  dynamical  systems  approach. 


The  Orr-Soinmerfcld  equation  is  discretized  using  a  finite  expansion  of  fi(i/)  (with 
(0,oo)  mapped  to  [— 1,+1]  using  an  algebraic  transformation)  in  a  Chebyshcv  series  re¬ 
ducing  (2.3)  to  an  algebraic  nonlinear  in  the  parameter  eigenvalue  problem  (see  Bridges 
&  Morris  [19S7])  for  this  type  of  reduction).  The  differential  eigenvalue  problem  (2.3) 
is  then  reduced  to  the  algebraic  eigenvalue  problem 


D4(A){fi}  =  [CoA4  +  Ci  A3  +  C2A2  +  C3A  +  C.,]{fi}  =  0.  (2.4) 


Note  that  in  the  fixed  frequency  spatial  stability  problem  (Bridges  &  Morris  [1987)) 
the  matrices  C0,...,C.i  are  complex  but  with  fixed  wavespecd  c  £  R  the  matrices 
Co, . . . , C.)  are  real.  The  eigenvalue  will  in  general  be  complex,  but  with  real  matrices 
the  numerical  computation  can  be  carried  out  with  greater  efficiency.  The  eigenvalue 
problem  (2.4)  is  solved  using  the  methods  of  Bridges  &  Morris  (1984)  for  nonlinear  in 
the  parameter  eigenvalue  problems. 

Given  (r,  R)  the  eigenvalue  problem  (2.4)  produces  a  spectrum  of  spatial  eigenval¬ 
ues.  The  idea  in  spatial  bifurcation  theory  is  to  look  for  bounded  solutions  of  (2.3);  that 
is,  there  arc  admissible  spatial  states  bifurcating  from  the  Blasius  boundary  layer  if  and 
only  if  there  exists  an  eigenvalue  A  of  (2.4)  with  Rc(A)  =  0.  Eigenvalues  with  Ro(A)  ^  0 
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arc  not  admissible  as  bifurcation  points  for  spatial  states  because  they  are  unbounded 
iis  *  — >  -|  oo  or  x  — *  —oo.  Note  that  we  are  not  concerned  with  eigenvalues  that  lie 
ofT  the  imaginary  axis  (whether  in  the  right  or  left  half  plane).  Ultimately  stability  is 
determined  by  checking  the  initial  value  problem  (reintroducing  time!).  If  Ro(A)  =  0 
and  lm(A)  =  0  there  is  (potentially)  a  bifurcation  to  a  new  equilibrium  state  (although 
this  is  not  expected  to  occur  for  the  Blasius  boundary  layer)  whereas  if  llc(A)  =  0  and 
Im(A)  0  there  is  a  (spatial)-IIopf  bifurcation  (assuming  the  usual  non-degeneracy 
conditions  on  Hopf  bifurcation)  from  the  equilibrium  state  to  a  spatially  periodic  state. 

Solving  |D4(A)|  =  0  and  Re(A)  =  0  results  in  the  well  known  neutral  curve;  for  the 
Blasius  boundary  layer  shown  in  Figure  2.1  (although  it  is  usually  plotted  in  (w,  R) 
or  (tv,  R)  space).  In  particular  there  is  at  most  one  pair  of  eigenvalues  on  the;  Im(A) 
axis.  Figure  2.3  shows  an  example  of  the  spectrum  of  the  eigenvalue  problem  (2.4) 
(there  is  also  a  (stable)  continuous  spectrum  of  (2.3)  (Grosch  &  Sahvcn  (1978])  that 
will  appear  as  discrete  in  the  finite-dimensional  approximation).  Accurate  calculations 
of  the*  ( c,R )  curve  with  the  associated  value  of  a  arc  given  in  Table  2.1.  Note  that 
there  is  a  finite  interval  in  wavespeed  c  G  (ci,C2)  in  which  the  neutral  curve  exists 
where  c\  ~  .22  and  C2  ~  .401.  We  call  c  admissible  if  c  G  (ci,c2).  It  is  clear  from 
Figure  2.2  that  if  c  G  (ci,C2)  and  fixed  and  R  is  increased,  a  spatial-) topf  bifurcation 
occurs  as  R  c tosses  the  neutral  curve.  There  is  a  continuum  of  Ilopf  bifurcation  points 
(varying  r)  and  this  will  have  consequences  with  regard  to  stability  (the,  possibility  of 
sideband  instability  will  have  to  be  considered)  but  nevertheless,  for  fixed  c  £  (cj,C2)  a 
classic  Hopf  bifurcation  takes  place  as  R  crosses  the  neutral  curve.  We  call  it  a  spatial 
Hopf  bifurcation  because  the  “frequency”  associated  with  the  bifurcation  is  in  fact  the 
slreamwise  wave  number. 

2.1  Primary  spatial  Hopf  bifurcation 

The  set  of  equations  (2.1)  can  be  written  as  an  evolution  in  x  in  the  following  way. 
Introduce  tin;  stream  function  with  u  =  ij>y ,  v  =  -ipx  and  the  vorlicily  £  =  -A i/>. 
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Then  the  vorticity  and  stream  function  equation  set  can  be  written  as 

n 

|U  =  L(c,7i)-«»  +  N($;rt) 


/V>\  /  V' 

,T,  _  I  V  I  'l£f  |  “V'r 

U  "  * 


,  N(d>, /?)  =  /?. 


W'»  + 


(2.0a) 


/  0  -1  0  0  \ 

Wa«)=  %  0  0  1  •  <2lil> 

V  o  ~nu"  -fa  mu-,-)} 

Taking  *1*  =  results  in  the  eigenvalue  problem  L(c,  r)'I»  =  A*l»  which  reduces 
to  the  Orr-Somtncrfcld  problem  (2.3)  for  the  stream  function  perturbation.  Suppose 
c  6  («: i ,  c-2 )  and  11  =  /?„  is  a  point  on  the  neutral  curve.  Then  ltc(A)  --  0,  lm(A)  =  a„ 
and  there  exists  an  eigenfunction 

L(c,/?)d>  =  ia0iK  (2.7) 

It  is  now  straightforward  to  apply  the  Ilopf  bifurcation  theorem  using  a  formal  centre- 
manifold  reduction  (Coullet-Spicgel  (1983]).  Scale  x  ax  so  that  the  wavenumber 
appears  in  c<|uation  (2.5):  o*Pr  =  L(c,  !?.)•$  +  N('I',  /?,).  Write  any  solution  of  (2.5)  as 

<I»(a:,;/)  =  A{x)<i'(y)  +  A(x)i>{y)  +  <fl(x,y) 

then  the  scaled  version  of  equation  (2.5)  is  transformed  to 

ilA 

aTl ~  =  fi(A,A,<i)  (2.8  a) 

=  (iM,¥)  (2.8a) 

a—  =  h{A,A,<i).  (2.8c) 


At  least  locally  (2.fc)  can  be  solved  for  'I'  as  a  function  of  A  and  A  (Coullet  fc  Spiegel 
[1983]).  Then  back  substitution  of  'P  into  (2.8a)  results  in  a  vcctorfield  on  C.  Applica¬ 
tion  of  normal  form  theory  (see  Guckcnhcimer  &  Holmes  [1983,  p.142])  results  in  the 
(formal)  normal  form  for  the  spatial-llopf  bifurcation, 


~  =  A-F(R-R0,a-aoM\2)- 
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The  reduction  of  (2.5)  to  (2.9)  is  ronual  simply  because  the  Blasius  solution  is  not  a 
truly  parallel  solution  of  the  Navicr-Slokes  equations  and  we  have  neglected  the  non¬ 
parallel  terms.  However  the  basic  idea  of  a  centre-manifold  reduction  for  a  spatial 
bifurcation  can  be  rigorously  justified  (Vandcrbauwlicdc  &  Iooss  [1990])  when  an  exact 
equilibrium  state  is  used  such  as  plane  Poiseuille  £  >w  (Bridges  [1991c]).  The  actual 
numerical  calculations  required  for  the  normal  form  (2.9)  will  be  discussed  in  Section 

3.2  and  aie  contained  in  Bridges  [1991b]  ns  a  special  case  of  the  3D  calculations. 

lixpaud  F  in  a  Taylor  series, 

F  =  F?t(R  -  R0)  +  F"(a  -  a0)  +  F“7|,l|2  +  •  •  • .  (2.10) 

The  imaginary  part  of  (2.10)  is  solved  for  (cv  -  a0)  and  back  substitution  into  the  real 
part  of  (2.10)  results  in  the  usual  pitchfork  bifurcation 

^  =  ag (11  -  Roy ),  a  =  g°n(R  ~  F„)  +  <7>2  -!-•••. 

Supposing  g'h  >  0,  a  supercritical  pitchfork  bifurcation  occurs  if  g„i  <  0  and  a  sub- 
critical  pitchfork  bifurcation  occurs  if  gni  >  0.  Stability  of  the  bifurcating  states  does 
not.  follow  directly  from  (2.9)  when  the  bifurcation  is  supercritical.  For  a  satisfactory 
stability  analysis  time  must  be  reintroduced  and  sideband  instability  considered. 

From  the  physical  point  of  view  the  normal  form  provides  two  crucial  pieces  of 
information:  (a)  the  direction  in  parameter  space  (Reynold’s  number  in  this  case)  in 
which  the  nonlinear  spatially  periodic  states  persist  ( R  <  Ra  or  It  >  R0)  and  (b) 
how  the  wavenumber  changes  along  a  branch  of  periodic  stales.  We  are  particularily 
inlet  ested  in  whether  the  wavenumber  goes  to  zero  along  the  branch  leading  to  spatial 
homoclinic  bifurcation. 

2.2  Secondary  bifurcations  and  spatial  Floquet  theory 

One  of  the  central  features  of  the  fixed-wavespeed  spatial-bifurcation  theory  is  that 
it  is  clear  how  more  complex  spatial  bifurcations  can  arise  (and  indeed  arc  expected). 
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In  Section  2.1  we  showed  that  for  for  fixed  c  €  (01,02),  a  spatial- 1 1  opf  bifurcation 
leads  to  a  branch  of  spatially  periodic  solutions  which  is  entirely  analogous  to  a  branch 
of  periodic  orbits  in  a  finite-dimensional  dynamical  system.  If  we  follow  a  branch  of 
periodic  solutions  then  we  can  expect  to  encounter  period-doubling  and/or  secondary 
bifurcation  to  a  quasi-pcriodic  flow.  The  idea  here  is  to  use  Floquet  theory  in  space; 
that  is,  to  study  the  spatial  Floquet  multipliers  along  the  branch  of  spatially  periodic 
solutions.  Let  (u,n,p)  (u  +  £,u  +  )/,j>+<j)  with  («,v,p)  a  periodic  solution  satisfying 
(2.1).  In  the  usual  way  substitute  (u  +  £,v  +  j/,p  +  g)  into  the  Navier-Stokes  equations 
and  liuentize  about  the  branch  of  periodic  solutions.  The  result  is  the  following  system 
with  periodic  coefiicients 


Ox  +  Oy 

(«  + ' V  -  C)  jg  +  «.£  +  (V,  +  ■ +  |2  -  Ia{  =  0 
(„  +  V  -  1  =  0- 


(2.11) 


This  set  can  be  simplified  by  introducing  a  perturbation  stream  function;  let  £  —  d<j>/dy 
and  »/  =  --#</»/ dx  then  the  second  and  third  of  (2.11)  can  be  combined  into 


AA0  -  /?(«  +  U-  c)^A +  R(vxt  -  uxy)°^ 

+  K(Vyy  +  Uyy  ~  v*y)  jr;  =  0. 


(2.12) 


Now  <l>(x,y)  is  not  necessarily  periodic  in  x  but  by  Floquet’s  theorem 


<^(x,y)  =  e7I^(x,t/)  7SC 


and  </>(*,//)  is  periodic  in  x  of  the  same  period  as  the  basic  state  (n,v,p).  Substitution 
of  <j>  =  c7*</)  into  (2.G)  results  in  the  following  eigenvalue  problem  for  the  Floquet 


exponent  7, 


L  •  4>  =  74L0^  +  7,’1Li<£  +  +  7L3<^  +  7L.i^  =  0 


(2.13) 
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where 


Lu  (j>  =  f> 

L,<£  =  4  jjt-R(u  +  U-c)j> 

L4  =  2A^  +  4^|  -  3/{(u  +  U-c)^-  Rv^ 

W  =  4Ag  -  /?(«  +  £/  -  c)  ^A*  +  20j  -  2i?„J^-  -I-  J?(VW  I-  «#y  -  „„)* 

W  =  AA<£  -  7?(u  +  u-  c)A^  -  /?.uA  +  R{Uyy  +  «,,,,  -  nr„) 

Not<‘  (hat  the  eigenvalue  problem  (2.13)  is  a  nonlinear  in  the  parameter  eigenvalue 
problem  for  the  Floquet  exponent  7.  It  is  reminiscent  of  the  classic  spatial  stability 
eigenvalue  ptoblem;  in  fact  it  is  the  spatial  form  of  the  secondary  “instability”  pioblem. 

The  above  theory  for  secondary  bifurcation  using  Floquet  theory  is  similar  to 
Herbert’s  [1983,1984]  theory  for  subharmonic  bifurcation  but  there  is  a  subtle  difference. 
The  eigenvalue  problem  (2.13)  is  the  spatial  secondary  “instability”  problem  whereas 
Herbert’s  theory  is  a  temporal  secondary  instability  theory.  In  Herbert’s  theory  the 
spatial  Floquet  exponent  is  treated  sis  fixed  (generally  so  that  the  spatial  multiplier 
lies  at  -J )  and  the  temporal  exponent  is  solved  for.  In  (2.13)  the  temporal  exponent  is 
absent  since  vie  arc  looking  for  secondary  steady  states;  that  is,  states  that  move  at  the 
given  spend  c  but  have  more  complex  spatial  structure. 

Given  a  Fourier-Chebyshev  representation  for  the  basic  spatially  periodic  slate, 
(a, a,/)),  the  eigenvalue  problem  (2.13)  can  be  discretized  by  expanding  in  a  finite 
Fom  icr-Ghebyshcv  series.  The  result  is  an  algebraic  nonlinear  in  the  parameter  eigen¬ 
value  pioblem 

(D..(7)] '  {*}  =  [CoY1  T  C1T3  +  C,72  -h  C37  +  C.,)  ■  {^}  =  0  (2.14) 

with  the  Floquet  exponents  obtained  as  roots  of  |D^(7)|  =  0.  Numerical  methods  for 
nonlinear  eigenvalue  problems  of  the  type  (2.14)  can  be  found  in  bridges  &  Mon  is 
[19S 1]  and  l’caristein  &  Goussis  [198S]  and  references  therein.  Given  the  spectrum  of 
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exponents,  the  Floquet  multipliers  are  given  by  exp(2;r7/c*).  Suppose  that  the  branch 
of  spatially  periodic  solutions  is  parametrized  by  a  parameter  6.  Then  two  interesting 
possibilities  as  c  is  varied  are  shown  in  Figure  2.3.  We  arc  not  concerned  with  where 
the  bulk  of  the  spectrum  of  die  eigenvalue  problem  lies  in  the  complex  plane  (although 
most  of  the  multipliers  will  probably  lie  in  the  interior  of  the  unit  circle);  any  multiplier 
not  on  the  unit  circle  is  inadmissablc  as  a  bounded  spatial  state.  Thercfoie  wc  vary 
e  until  the  multiplier  lies  on  the  unit  circle.  If  the  multiplier  passes  through  —1  we 
expect  a.  wavelength  doubling  bifurcation  and  if  complex-conjugate  multipliers  pass 
through  the  unit  circle  at  points  other  than  ±1  we  expect  a  bifurcation  to  a  spatially 
qunsipet  iodic  state. 

Numerical  calculations  arc  necessary  to  obtain  precisely  where  secondary  bifurca¬ 
tions  occur.  However,  spatially  quasi- periodic  states  are  to  be  expected.  One  way  to 
show  this  is  to  introduce  a  singularity  (in  the  neutral  curve)  which  results  in  a  larger 
dimensional  centre-manifold,  then  look  for  complex  dynamics  in  the  unfolding.  In  fact, 
in  Section  5  it  is  shown  that  the  compliant  wall  problem  has  a  .singularity  or  this  form 
fiom  which  we  can  show  the  existence  of  secondary  branches  of  spatially  quasi-periodic 
stales  (and  possibly  spatially  quasi-periodic  states  with  three  frequencies!). 

The  wavelength  doubled  solution  can  again  double  its  wavelength  with  a  possible 
cascade  to  spatially  chaotic  states  (not  turbulent  though  since  wc  arc  restricted  to  two- 
dimensions  but  the  3D  problem  is  considered  in  Section  3).  An  interesting  theory  for 
wavelength  doubling  cascades  can  be  constructed  as  follows.  Note  that  from  Figure  2.1 
that  for  fixed  c  €  (ci,C2),  there  is  only  a  finite  band  in  R  in  which  ail  unstable  region 
occurs.  From  numerical  calculations  of  Herbert  [1975]  we  expect  the  nonlinear  branch  of 
periodic  solutions  lo  form  a  global  loop  as  shown  in  Figure  2.4.  The  global  wavelength 
doubling  structure  of  loops  can  be  modelled  by  a  one-dimensional  two-parameter  map. 
For  example  consider 

*„+i  =  *H  (1  -  (R  -  Ro)2  +  7  -  *?,)•  (2.15) 

The  fixed  points  of  the  map  (2.15)  are  as  shown  in  Figure  2.4(a).  Iterating  the  map 
while  increasing  7  (corresponds  to  decreasing  the  wavespeed  c)  results  in  secondary, 
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tertiary,  etc  loops  of  period  doubled  points  (Bridges  (1091ej).  Figure  2.5  shows  a  finite 
cascade  of  wavelength  doubling.  If  7  is  increased  further  the  cascade  will  become  infinite 
resulting  in  chaos  in  a  thin  subinterval  of  R.  It  is  very  likely  that  this  is  the  structure 
of  wavelength  doublings  in  shear  Hows.  Secondary  period-doubling  bifurcations  from 
the  subcritical  loop  in  Figuic  2.4(b)  can  be  modelled  in  a  similar  fashion  using  the 
one-dimensional  map 

H  I  =  a.'n(l  -  (R  -  R0)2  +  y-(ft  +  2 m(R  -  R0))r2  -  z*u)  (2.16) 

where  /J,»n  €  R,  0  <  m  <  1,  7  >  0  and  ft  <  —  >/4m27.  The  period  doubling  structure 
for  the  map  (2.1G)  is  shown  in  Figure  2.6  with  a  finite  cascade  in  Figure  2.6(a)  and 
an  infinite  cascade  leading  to  a  region  of  chaos  in  Figure  2.6(b).  Note  that  although 
the  subcritical  and  supercritical  loops  in  Figure  2.4(a)  and  (b)  differ  significantly,  the 
secondary  structure  and  cascade  structure  is  about  the  same  in  both  cases. 

Note  that  the  above  theory  is  restricted  to  two  spatial  dimensions.  If  indeed  spatial 
chaos  occurs  it  will  not  lie  turbulence.  The  role  of  three-dimensionality  is  considered  in 
Section  3.  On  the  other  hand,  three-dimensionality  does  not  affect  the  cascade  theory 
shown  in  Figure  2.5,  it  simply  results  in  non-trivial  spanwisc  structure  along  with  the 
streamwise  cascade. 


2.3  Spatial  evolution  of  the  primitive  variables 


As  an  alternate  to  the  evolution  equation  (2.5)  where  the  stream  function  vort ic¬ 
ily  variables  arc  used,  an  evolution  c  .nation  for  the  primitive  variables  can 

be  constructed  although  it.  has  a  nonstandard  form  (involving  an  evolution  equation 
and  a  constraint).  The  idea  is  to  evolve  (in  x)  the  pressure  and  vertical  velocity  and 
determine  the  streamwise  velocity  using  a  constraint  (a  differential  equation  without 
^-derivatives).  First,  construct  a  Poisson  equation  for  the  pressure  by  taking  the  diver¬ 
gence.  of  the  momentum  equations  in  (2,1), 


Ap  =  -2 


f  dv  du  dv  dv  dv\ 
\  v  dx  ^  dy  dx  ^  dydy) 


(2-17) 
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Now  lei. 

:)• 

<1  / 

Thou  tlio  Poisson  equation  (2.17)  and  the  vertical  momentum  equation  can  together 
be  written  as  an  evolution  equation  in 

r\ 

=  L(c,  R)-i>  +  N($,«;  R)  (2.19) 


where 


( 

0 

1 

0 

0\ 

(  0  \ 

Me,  R)  — 

— 

a1 

7h,r 

0 

R{U  -  c) 

0 

rM. 

P 

« 

1 

and  N  = 

f?(uVM-  .»g) 

0 

\ 

0 

"2  Uy 

1)7 

V-2(g)2-2  VfJ 

Note  however  that  u  appears  in  the  nonlinear  term  but  is  not  a  component  of  $. 
However,  the  streamwise  momentum  equation  can  be  written  as  a  differential  equation 
in  t)  alone, 

<1  +  UyV  -  ( V  -  c)v,j  +  i(Vy  -  it„„)  -  uvy  +  vuy  =  0,  (2.21) 

and  at  each  value  of  *,  it  is  obtained  from  (2.21)  for  use  in  (2.20). 

The  evolution  equation  for  («,  u,p)  is  non-standard  in  that  it  involves  evolution  of 
( t»(/<)  with  a  constraint  to  determine  the  streamwise  velocity.  Nevertheless  it  is  a  useful 
form  of  the  equations;  in  particular,  the  framework  (2.19)-(2.21)  is  easily  extended  to 
the  3-dimensional  Navior-Stokes  equations  (this  is  carried  out  in  Section  3.5)  and  the 
usual  ceuti e-manifold  and  bifurcation  theory  is  still  applicable. 


(2.18) 
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Table  2.1 


Coordinates  in  the  (c,it)  plane  and  wavenumber 
along  the  neutral  curve  for  the  Blasius  boundary  layer 

Upper  branch  Lower  branch 


R 

c 

Q'o 

R 

c 

Q'o 

4000 

.2956 

.2689 

4000 

.22875 

.0918 

3500 

.3028 

.2777 

3500 

.2368 

.0967 

3000 

.31065 

.2881 

3000 

.2461 

.1028 

2600 

.3185 

.2976 

2560 

.2552 

.1097 

2500 

.3206 

.3001 

2500 

.2568 

.1106 

2400 

.3229 

.3029 

2360 

.2608 

.1135 

2200 

.3276 

.3086 

2160 

.26655 

.1178 

2000 

.3329 

.3148 

2000 

.27105 

.12175 

1800 

.3387 

.3215 

1960 

.2726 

.1229 

1600 

.3451 

.3286 

1760 

.2798 

.1290 

1500 

.3489 

.3322 

1560 

.2881 

.1364 

1400 

.3528 

.3362 

1500 

.29085 

.1389 

1200 

.3616 

.3440 

1360 

.2981 

.1456 

1000 

.3718 

.3515 

1260 

.3034 

.1513 

900 

.3777 

.3545 

1160 

.3105 

.1578 

800 

.3843 

.3562 

1060 

.3178 

.1655 

750 

.3877 

.3562 

1000 

.3223 

.17095 

730 

.38915 

.3559 

960 

.3259 

.1749 

660 

.3941. 

.3532 

900 

.3313 

.1814 

650 

.3947 

.3524 

850 

.3363 

.1S79 

600 

.3981 

.3469 

800 

.34205 

.1950 

580 

.3992 

.3432 

750 

.3484 

.2035 

575 

.3993 

.3421 

700 

.3550 

.2135 

560 

.4001 

.3379 

650 

.3630 

.2259 

550 

.4004 

.3344 

600 

.37205 

.2419 

540 

.4002 

.3293 

580 

.37615 

.2499 

530 

.4002 

.3241 

560 

.3812 

.2597 

525 

.3995 

.3184 

550 

.38365 

.2656 

520 

.3951 

.2974 

540 

.3804 

.2725 

530 

.3987 

.2815 

525 

.3945 

.2S57 

n 


Figure  2.4  Global  loop  structure  for  fixed  c  of  spatially  periodic  states  bifurcating 
from  the  Blasius  boundary  layer:  (a)  supercritical  loop  and  (b)  subcritical  loop. 


Jl 


Figure  2.5  Global  period-doubling  loops  with  a  finite  cascade  in  the  map  (2.15):  (a) 
7  <  1  resulting  in  an  absence  of  period-doubling  and  (b)  7  >  1  (in  particular  7  =  1.30) 
resulting  in  three  period-doubing  bifurcations. 


Figure  2.6  Finite  and  infinite  period-doubling  cascade  in  the  map  (2.16)  where  the 
primary  loop  is  subcritical  with  in  =  5,  /?  =  -y/j-  .2  and  (a)  7  =  .21  and  (b)  7  =  .24. 
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3.  Spatial  bifurcations  in  three  dimensions 


In  this  section  we  continue  to  treat  the  steady-state  (in  a  moving  frame)  Navier- 
Stokes  equations  as  an  evolution  equation  in  the  streamwisc  coordinate  but  with  the 
addition  of  spamvisc  variation  (three  dimensionality).  The  basic  equilibrium  state  is 
taken  to  be  the  (parallel)  Blasius  boundary  layer  although  the  theory  is  generally 
applicable  to  any  two-dimensional  parallel  equilibrium  state.  Wi  ‘.h  the  shift  a;  t--»  a:  —  c t 
the  generalization  of  the  equation  set  (2.1)  to  three  dimensions  is 


du 

dx 


dv  dw 
dy  *  i h 


=  0 


...  .du  dp  1  .  du  du  du 

(if  -  «)S  +  u,v  +  m  +  ..jj  + » gj  +  =  » 

<»-«>£  +«S->A.  +  .*+,^  +  »^-0 

dx  dy  11  dx  cy  dz 

...  . dw  dp  1  .  dw  dw  dm 

w~c)0i  +fe-7iAn,+"fer't"’8?+"’or  =  0 


(3.1) 


The  addition  of  z  (the  spamvisc  coordinate)  introduces  nontrivial  symmetry  that 
will  be  fundamental  to  the  analysis  of  three-dimensional  states;  in  particular,  a  trans¬ 
lation  and  reflection  in  z.  The  translation  and  reflection  in  z  follow  from  the  fact  that 
the  basic  slate  is  two-dimensional.  The  --reflection  generates  the  group  Z$  “  (n)  with 
action  on  functions  given  by  n-  f(x,y,z)  =  f(x,y,-z).  The  set  (3.1)  is  ZJ-equivariant 
and  {u(x,y>-z),v(x,y,-z),-w(x,y,-z)tp(x>y,-z))  is  a  solution  of  (3.1)  whenever 
(u(x,y,z),u(x,y,z),w{x,y,z),p(x,y}z))  is  a  solution;  that  is.  there  is  a  Zj -orbit  of  so¬ 
lutions.  Note  that  the  existence  of  a  ZJ-orbit  of  solutions  is  true  regardless  of  solution 
type;  even  chaotic  trajectories  have  a  Z^’ -orbit  which  has  important  consequences  with 
regard  to  division  and  multiplicity  of  attracting  sets. 

The  translation  invariance  in  z  of  the  set  (3.1)  results  in  a  group  orbit  or  solutions 
as  well  (an  arbitrary  translate  in  2  of  a  solution  is  also  a  solution).  However  if  (u,  v,  w,  p) 
is  taken  to  be  periodic  in  z  (an  assumption;  more  complex  spamvisc  spatial  structure  is 
possible  and  this  is  considered  in  Section  3.4)  then  the  translation  group  is  reduced  to 
the  compact  group  SO(2).  Combining  SO(2)  with  Z£  results  in  an  0(2)-equivariance 
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of  the  equation  set  (3.1).  In  our  subsequent  analysis  the  0(2)-equivariance  of  the  set 
(3.1 )  is  I  lie  basic  organizing  feature  of  the  three-dimensional  states. 

In  three  dimensions  the  basic  spatial- Hopf  bifurcation  introduced  in  Section  2.1 
persists  but  the  0(2)-cquivariance  results  in  a  higher  dimensional  (spatial)  centre- 
manifold,  a  higher  multiplicity  of  spatially  periodic  states  and  the  potential  for  more 
complex  “dynamics”  (i.c.  more  complex  spatial  bifurcations).  For  the  spatial  Hopf  bi¬ 
furcation  with  0(2)  symmetry  we  adapt  the  well-developed  theory  of  ()(2)-cquivarinut 
Hopf  bifurcation  (Golubitsky  &  Roberts  [1987],  Golubitsky,  Stewart  &  Schaeffer  [1988]). 
The  0(2)-equivariant  spatial-Hopf  bifurcation  is  a  primary  bifurcation  to  3D  states  and 
is  lieated  in  Section  3.1. 

In  Section  3.3  we  introduce  a  “spatial”  secondary  instability  theory  where  a  pri¬ 
mary  spatially  periodic  two-dimensional  state  bifurcates  to  a  three-dimensional  state 
at.  finite  amplitude  in  a  steady  frame.  This  is  to  be  compared  with  the  temporal  sec¬ 
ondary  instability  theory  due  to  Orszag  &  Patera  [1983]  and  Herbert  [1983,1984].  The 
2D  ->  3D  spatial  secondary  “instability”  theory  is  similar  to  the  theory  intioduced  in 
Section  2.2  but  includes  nontrivial  (periodic)  spamvise  variation.  Linearization  about 
the  2D  stale  results  in  a  system  with  periodic  (in  x)  coefficients  to  which  spatial  Flo- 
quet  theory  is  applied.  The  system  will  differ  from  the  system  (2.11).  The  spalial 
Floquet  exponents  now  depend  on  two  parameters:  the  parametrized  branch  and  fi 
(the  spamvise  wavenumber  of  the  perturbation).  Our  approach  differs  from  the  tem¬ 
poral  theories  of  Orszag  &  Patera  and  Herbert  in  that  the  temporal  exponent  is  set 
to  zero;  that  is,  we  look  for  bifurcations  to  bounded  steady-states  with  moie  complex 
spatial  structure  in  x  (wavelength  doubled,  tripled,  etc.  or  quasi-periodic.)  but  periodic 
in 

It  is  clear  that  if  a  primary  state  exists  that  is  periodic  in  both  the  stroamwisc 
direction  and  the  spamvise  direction  (as  obtained  in  Section  3.1  or  as  a  secondary 
bifurcation  as  in  Section  3.3),  it  is  possible  to  have  secondary  bifurcations  in  both  the 
stieamwise  and  spamvise  directions.  In  particular,  it  is  possible  to  have  a  bifurcation  to 
states  with  more  complex  spanwise  spatial  structure.  The  theory  for  spatial  bifurcation 
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in  (a:,  z)  will  involve  spatial  Floquet  theory  in  two  space  dimensions  and  is  discussed  in 
Seel  ion  3.4. 


3.1  0(2)-equivariant  spatial  Ilopf  bifurcation 


To  determine  bifurcation  points  from  the  equilibrium  state,  the  set  (3.1)  is  lin¬ 
earized  resulting  in 

Ou  Ov  Out 

=0 

which  can  l>c  reduced  to  the  two  decoupled  systems 

Av  -  R{U  -  c)^  -  I&-  =  0 

dx  dy  (3.3) 

Ap  +  2t/j,~  =  0 
Ux 

ami 

An  -  R(U  -  c)^  =  R^-  +  RUyV 

'Ox  Ox  (3.4) 

a  li/rr  \°w  r>VV 

Am  R{U-c)dx  =Rgz- 

Taking  «(.t,j/,z)  =  eAl[ai(y)cos/iz  +  Ui(y)  sin  Pz]  the  secondary  eigenvalue  problem 
(3.4)  reduces  to 

^  +  [A2  -  P2  -  A R(U  -  =  0  i  *  1,2  (3.5) 

with  appropriate  boundary  conditions.  It  is  easy  to  show  that  if  <:  G  11  (and  R  finite) 
tluat  cocry  member  of  the  point  spectrum  of  (3.5)  is  real  and  non-zero.  Potential 
bifurcation  points  are  therefore  obtained  from  the  eigenvalue  problem  (3.3).  Let 


teO-to[tesi)“"+(sso 


sin/?z  , 


tl»  n  (;&  -I-  (A2  -  P2))Vj  +  2\UyVj  =  0  (j  =  1,2)  and 

-I-  (A2  -  P2)  1  Vi  -  A R(U  -  c)  ~  +  (A2  -  p2)  Vi  +  \RUy,jVi  =  0  (3.7) 
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which  is  a  modified  (real)  form  of  the  three-dimensional  Orr-Sommerfeld  equation. 
Given  (e, It,/))  (all  real)  there  are  two  linearly  independent  eigenfunctions  for  each 
eigenvalue  A  of  (3.7),  hence  each  eigenvalue  of  (3.7)  has  (gencrically)  geometric  and 
algebraic,  multiplicity  two.  Consequently  if  there  exists  an  eigenvalue  A  of  (3.7)  with 
Re(A)  —  0  and  Iin(A)  ^  0  then  it  is  (gencrically)  double  and  with  its  complex  conjugate, 
the  associated  4  eigenfunctions  form  a  basis  for  a  four-dimensional  (spatial)  centre- 
subspace. 

Hopf  bifurcation  points  arc  found  by  fixing  (c,/?)  and  increasing  11  until  there 
exists  an  eigenvalue  of  (3.7)  with  llc(A)  =  0  and  Im(A)  ^  0,  or  for  fixed  ft  the  neutral 
curve  in  the  (e,  R)  plane  can  be  obtained.  When  ft  =  0  the  neutral  curve  is  as  shown 
in  Figure  2.1  and  when  ft  ^  0  Squire’s  theorem  (Drazin  &  Iteid  (1981,  p.  155])  can  be 
adapted  to  determine  the  ft  ^  0  neutral  curve:  suppose  ft  =  0  and  A  =  m  (a  €  R) 
and  let  (c,  /?„)  be  a  point  on  the  neutral  curve  with  wavenumber  n„.  Then  for  ft  0 
Squire’s  theorem  states  that  for  each  given  c  the  neutral  point  (for  ft  0)  is  shifted 
(positively)  to  Rp  =  Ron0/otp  where  op  =  \/al  —  ft2  (assuming  \ft\  <  |o0|)-  Given  the 
neutral  curve  for  ft  =  0  it  is  therefore  straightforward  to  construct  the  ft  ^  0  neutral 
surface  and  curves  for  various  values  of  ft  are  shown  in  Figure  3.1. 

Nole  that  for  each  admissable  ft  there  exists  a  point  in  (c,  R)  space  at  which  pure 
imaginary  eigenvalues  for  the  2D  and  3D  states  exist  simultaneously  (where  the  ft  =  0 
and  ft  •-]-  0  curves  intersect).  These  points  are  codimension  2  bifurcation  points  and  since 
the  eigenvalues  arc  gencrically  nonresonant,  they  correspond  to  points  of  bifurcation  to 
spatially  quasi-periodic  states  and  arc  analyzed  in  Section  4.  For  particular  values  of  ft 
the  eigenvalues  of  the  codimension  2  points  will  be  resonant  (this  will  be  of  codimension 
3).  For  example  if  at  ft  =  ft0  the  2D  and  3D  eigenvalues  lie  at  a0  and  £cv„  the 
lesonance  is  a  spatial  version  of  the  Craik  resonant  triad  (Craik  [1971,1985]).  In  Section 
4  these  resonances  (o0  and  or0/n  n  =  2,3,4)  are  considered  from  the  spatial  point  of 
view;  that  is,  they  are  codimension  three  organizing  centers  for  more  complex  spatial 
structuie.  In  this  section  it  is  assumed  that  (c,ft)  take  generic  (admissable)  values  and 
the  codimeiision  1  bifurcations  associated  with  variation  of  11  are  treated. 
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Wo  now  proceed  to  compute  the  bifurcation  from  the  neutral  curve  with  /?  ^  0. 
Choose  admissablc  values  of  /?  and  c  and  suppose  R  =  R0  is  the  neutral  point  with 
eigenvalue  A  =  ia„  with  a0  6  R.  Scale  a:  t-*  a  a  so  that  a  appears  as  a  coefficient  in  the 
nonlinear  set  of  equations  (3.1).  Reverting  to  complex  coordinates  the  solution  of  the 
linear  equation  (3.2)  at  the  neutral  point  is  given  by 

«i(*,2 /.2)\  [  /  Mv)\ 

Vi (x,y,z)  =2Rc  (Ae'l'+W  +  Z?cl(*-^*>)  I  0,  (j/) 

vi(x,y,z)  J  [  \Mv)J. 

and 

w>i(*,  yfz)  =  2Re  [( Aci<r+'}**>  -  (t/)] 

whore  A,D  €  C  arc  complex  amplitudes.  Let  A  ^  —  (a2  -I-  /?2)  then  the  complex 
functions  (»i,V2,«>3,j>.t)  satisfy 

U{a„0)  •  vi  =f  A2S,  +  ia0RUyyVi  -  ia0R{U  -  c)Af>,  =  0  (3.9) 


m0  dvi  t  f)2R  T_,,rr  „  , 

« J  =  -?7  en  ~T~  +  -  ?  ;  ^,-L;  (t/jUl) 
al  +  132  Uy  a2  +  /?2 

.  __  W  di>  1  OCoftR  T-i/TT  ~  N 

WI  Ql+U^2  {UyVl) 


(3.10) 


whore  L2  *=  A  -  ia0R(U  -  c). 


The  idea  is  to  apply  the  ecu  tic-manifold  theorem  to  the  spatial  bifurcation  of 
periodic  states.  Rigorous  application  of  the  centre-manifold  theory  is  not  possible  in 
this  case  due  to  the  fact  that  the  Blasius  solution  does  not  satisfy  the  Navier-Stokcs 
equations  and  the  further  neglect  of  non-parallel  terms.  We  can  however  give  a  formal 
construction  of  the  centre-manifold  using  the  theory  of  Coullet  &  Spiegel  [1983].  There 
are  two  steps  in  the  reduction  to  normal  form.  Suppose  (3.1)  has  been  recast  iis  an 
evolution  equation  (this  is  carried  out  for  the  primitive  variables  in  Section  3.5), 


L(c,i?)  •  <I>  +  N(4>,m,  ;7l). 


(3.11) 
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For  adiuissablc  c  and  R  —  R„  the  linear  operator  L  satisfies 


L(c,  R)  •  !>(»/,  r)  =  io0^(t/,  z)  a0  €  R 

where  z)  =  (Ae,l,!  +  Bc_,/K')'i'i(j/)  with  A,#  G  C.  Suppose  that  A  =  dkt’cr,,  are 
the  only  eigenvalues  of  L(c,  R0)  on  the  imaginary  axis.  Then  <I>  can  l)C  expressed  as 

‘I >(*.».*)  =  {A(x)cil,:  +  D(x)c-^yh(y) 

+  (A&jc-#*  +  +  *(x,y,  z)  (3.12) 


where  'I'  consists  of  all  modes  not  in  the  center  subspace.  The  first  step  in  the  centre- 
manifold  reduction  is  to  substitute  (3.12)  into  (3.11)  (at  R  =  R0 )  resulting  in 


-£:  =  Ma,b^,b,*) 

~  =  Ma,b,a,b,<z) 


(3.13) 


(^  =  f3(AB,A&*)  (3.14) 

with  additional  equations  for  A  and  B.  All  eigenvalues  of  (//3(0,  (),(),  0,0)  are  ofr  the 
imaginary  axis  and  presumably  (3.14)  can  be  solved  for  $  :ts  a  function  of  AfB,  A>B. 
Hack  substitution  of  4'  into  (3.13)  results  in  a  vectorfield  on  C2.  The.  reduced  vcctorfiold 
will  be  an  0(2)-cquivariant  vectorfield  on  C2  with  ±»a0  (double)  eigenvalues  at  the 
linearization.  Normal  form  theory  is  then  applied  to  transform  (3.13)  to  an  0(2)  x  S1- 
equivarianl  vectorfield  on  C2; 


~  =  AM«-a0,R-RoM\2M) 

^  =  BMa-a0,R-R0,\A\2AB\2) 

but  the  0(2)  x  S1  symmetry  requires  /^(v, |A|2,  |i?|2)  =  /)(•,•, |#|2, |/1|2). 


(3.15) 


First  we  will  analyze  the  normal  form  equations  (3.15)  using  the  theory  of  Gol- 
ubilsky,  Stewart  &  Schaeffer  (19SS,  Chapter  XVII]  and  Golubitsky  &  Roberts  [1987] 
keeping  in  mind  that  the  “frequency”  is  in  fact  the  wavenumber  a.  Then  details  of  the 
construction  of  (3.15)  are  given. 
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Lcl.  N  =  |/1|2  +  |R|2,  A  =  S2  and  8  =  |i?|2  -  |/1|2  then  following  Golubitsky  & 
Roberts  [1985,  prop.  2.1)  the  set  (3.15)  can  be  transformed  to 

^4  =  (p  +  iq  +  (r  +  ia)6)A 

(3.16) 

D  =  (p  +  iq  —  (r  +  is)8)D 

where  p.  q,  r  and  s  are  0(2)  X  S1  invariant  functions;  that  is,  they  arc  functions  of  N,  A 
and  parameters.  Writing  A  =  ae'^'1  and  D  =  be"!’2  the  set  (3.1G)  can  be  split  into 
amplitude/ phase  equations 

a  =  (p  +  r8)<  t  'j 

>  amplitude  equations  (3.17) 

b  =  (p-v8)b  J 

rj>l  ={q  +  s6)  1 

>  phase  equations.  (3.18) 

ii>2  =  (q-  *8)  J 

The  idea  is  to  solve  the  phase  equations  for  a  -  a0.  Then  substitution  of  the  expression 
for  tv  -  n0  into  (3.17)  results  in  a  function  of  71  —  R0  and  (a,  6)  alone  (when  («,/?)  are 
fixed): 

where 

y(a,b,R-  R0)  =  p{N,A,R-  R.„)  j  +  r(W,A,Jl-  R„)$  (  “fcJ  .  (3.19) 

Generieally  there  are  two  types  of  solutions  of  the  normal  form:  (a)  oblique  waves  with 
«  0  and  6  =  0  which  correspond  to  waves  with  a  wavefront  at  some  angle  to  the 

streamwise  direction  and  (b)  “standing  waves”  with  a  =  6.  These  correspond  to  waves 
that  travel  in  the  streamwise  direction  but  arc  periodic  in  the  spanwiso  direction.  There 
is  also  I  he  possibility  (with  an  additional  parameter)  for  the  two  classes  of  waves  to 
interact  producing  quasi-periodic  waves.  Further  details  on  the  symmetries  and  more 
complete  analysis  of  the  normal  form  (3.1G)-(3.19)  can  be  found  in  Chapter  XVII  of 
Golubitsky,  Stewart  &  Schaeffer  (19SS). 
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The  coefficients  in  the  normal  form  (3.12)  are  obtained  in  the  following  way.  The 
eigenfunction  4»i(»/)  in  (3.12)  is  easily  constructed  using  the  expressions  in  (3.8)-(3.10). 
The  equation  (3.14)  is  solved  only  to  sufficient,  order  (in  A  and  13)  to  obtain  the  generic 
normal  form.  This  is  done  by  expanding  $  as  a  quadratic  polynomial  in  j4,R,7T  and 
11.  Dropping  terms  that  don’t  appear  in  the  normal  form,  is  constructed  from 
(n*2 ,  a 2 ,  ie*2 )  where 


v2{x,yyz)  =  2Re  [(A*eW(,-w*)  +  B2e2i{t~fi,))  hi(y)  +  2ABclirvu{y)\ 
+  (AD'cr,lU  +A*Bc-Wl)  v23{y) 


(3.20) 


where 


I  _  A  nd  ( (Hi i  dih  A  d2vi\ 

--L,(-o .,-/J)-t»3I  -  ) 

^-L.t(2ao,0)  •  W22  =  4al~(u]  -  v\)  -  2 ia„  +  4a2^  ft jt>i 

~L, (0,2/3)  •  t)23  =  -8/?*^(M2  +  |«>.|2)  -  2 ip  (~  +  4/32)  (fljtfii  -  MJ).  j 

(3.21) 

For  the  strenmwi.se  and  spamvise  velocity  we  find 


«2(*,»,*)  =  2Re  [(J4V'(j+"i>  +  B 2e2,'(J-/J‘‘))u2i(y)  +  2>llie“fl5«22(?/)] 

+  (^*e2f"s  +  A'Be-W’YiM  +  (M|2  +  |73|2)M.V)  (3.22) 
»t>2(*,tf,2)  =  2Rc  [(yl2c2,'(,+/,-'J  -  Z?2e2f<*-^*)),&2i(lf)] 

+  (AB*c2ill:  -  A*Bc-2il,*)tV23(y)  +  (\A\2  -  \B\2)xv2i(y)  (3.23) 
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where 


— =  2»<v0«2i  (.v)  +  2iftw2i{y), 

f^7  -  *10I)tt2i (a)  =  tlUfVis  -I-  R~{ui<ii  +  2i/?/I(«ii6|  - 

^«2t(.v)  =  -i-  njOi) 


^«»2i(»/)  =  R(w\vl  4-  iiiJOj).  j 

,!/  (3.24) 

Substitution  of  (3.20)-3.24)  into  the  vcctorfield  (3.13)  and  subsequent  normalization 
results  in 


~  =  4  (/»3|(ii- no)  +  /,a,(o  -  (.„)  +  /tjjMI2  +  MZ?|2  ■  •  •] 
•7  =  -4  [/iai(/i  —  Ro)  T  /i3'i(a  —  a0)  +  tarMI2  +  taal^l2  T  •  •  ■] 


(3.25) 


where  lisj  -  [<t>*Jl3j{ll)\  and  </'•  is  the  adjoint  eigenfunction  of  the  Orr-Sommerfold 
equation  and 


l>3t(y)  =  ^-AAfq, 

r  ,  »  2m0  ;  .  2  .(/iti  nT,  ,dv l  .  .  .  , 

/’.,(»)  =  +  2Cr,(-¥  +•«.«,) 

+  tfi/yf'i  -I-  (17  -  c)At>i-  +  2 ia0(U  -  c)- 

hz{y)  =  -I-  («2  +/?2)  kn(y)  + -^jkuiy) 

hn(y)=  -j~2  +  (ao  T  /^2)  hi[y)  +  j-k22{y). 


(3.26) 


The  operator  A  is  the  reduced  Laplacian,  A  '=f  ^  —  (a2  4-  /P)  and  the  functions  t,-y 
are  defined  in 


/-•ii  =  (— iovtj  -  iPw\)t>2\  -  (»«„« 21  +  ifiw2\)vi 
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-  if),*(2itte,«2i  +  2f/?t02i ) 

k 21  -=  (-lOfoiti  +  »/?«> i)f>23  +  2f>2i{ifhb\  -  icf0Uj)  -  2ia0vfu22 


(3.27a) 


+  vi(if)w-u  ~  ia0uu)  -  vi(iar0u33  +  !/?t?i23)  (3.276) 

fci2  =  — 2(ia0«i  +  i(hvi)(hx0U2  \  +  f/3««'2i)  -  2(a*  +  /?'2)u,  t)2i 

+  2(-ia0uJ  -  iftib\)(ia0ii2i  *1*  */iit>2i )  (3.27c) 

A’22  =  ~2(»tVo«i  +  »7?U>| )(far0«2 1  -  »/?»»?»)  -  2(ia0tti  -  i(irui)(ia„V23  +  t’/b&zs) 

- 4taB«22(*rto«i  -  t/hi'i)  -  2(cv'g  +  02)(vtv-23  +  25(t>,2)  (3.27. i) 


Tlic  veclorfield  in  (3.25)  can  be  recast  as 

■J7  =  -4  [/*3l(7?.  —  /?o)  +  632(0  —  Oo)  +  |(/t33  +  lhi)N  +  5(634  —  633)5  -)■•  ••] 


which  is  of  the  form  (3.1G)  with 

V  +  »7  =  631(7?.  -  7?o)  +  632(0-  -  o-o)  +  5(633  +  634)//  -I- 
r  4-  is  =  5(631  —  633)  H-  •  •  • . 


(3.2S) 


It  is  now  straightforward  to  apply  the  theory  for  0(2)  x  S1  cqui variant  normal  forms 
to  (3.28)  given  the  complex  numbers  Ihj  j  =  1,...,4.  Numerical  evaluation  of  the 
cocllicicuts  in  (3.2S)  is  considered  is  Section  3.2. 


It  is  important  to  note  that  the  above  bifurcations  are  spatial  bifurcations  and 
the  stability  assignments  given  in  Golubitsky,  Stewart  &  Schaeffer  [J  OSS,  Cliapl.  XVII] 
arc  not  applicable.  To  determine  the  stability  of  the  two  classes  of  waves  (oblique  and 
standing)  time  will  have  to  be  reintroduced  and  the  possibility  of  sideband  instability 
consideied.  This  is  a  very  interesting  problem  that  we  will  treat  in  detail  elsewhere. 


3.2  Compulation  of  the  coefficients  for  spatial  Hopf  bifurcation 

I11  this  Section  details  of  the  numerical  evaluation  of  the  coefficients  in  the  0(2)xS* 
cqui  variant  normal  form  obtained  in  Section  3.1  are  given  for  the  case  where  the  equilib¬ 
rium  state  is  the  Blasius  boundary  layer.  Numerical  evaluation  is  essential  because  the 


33 


Dlasius  solution  is  given  implicitly  as  the  solution  of  a  differential  equation.  The  basic 
problem  is  lo  construct  the  functions  lisj{y)  ( j  =  1,...,4)  in  equation  (3.26)  which  in 
turn  depend  on  complex  functions  that  satisfy  the  homogeneous  or  inhomogeneous  Orr- 
Sommcifcld  (or  related)  equation.  The  Dlasius  equation,  the  Orr-Sommcrfcld  equation 
and  related  equations  are  all  posed  on  the  semi-infinite  interval  y  6  (0,  oo).  However  nu¬ 
merical  solution  of  the  Orr-Sommcrfcld  equation  on  semi-infinite  domains  is  treated  in 
Bridges  &  Morris  [19S4b,1987]  and  we  use  their  basic  algorithm  here.  The  semi-infinite 
domain  is  mapped  to  [-1,+lj  using  the  algebraic  transformation 

y  »€(0,oo).  (3.29) 

y  +  * 

Then  i->  nt(Y)jf  where  m(Y)  =  (1  — 1'")2/4  and  all  functions  of  y  are  conshlered  as 
functions  of  Y  and  expanded  in  finite  series  of  Ghebyshev  polynomials.  As  an  example 
we  consider  the  construction  of  a  finite  Ghebyshev  series  expansion  for  the  complex 
function  v-n(y)  (in  equation  (3.19))  which  satisfies  an  inhomogeneous  (modified)  Orr- 
Sommcifcld  equation.  Mapping  y  *-*  Y  the  governing  equation  for  Dai  is 


+  /l(n»4+B(K)«n  = 


(3.30) 


where 


/l(y)  =  -2  (al  +  p)-iann0(U(Y)-c) 

H(Y)  =  («*  +  P)2  +  +  P2){U  -  c)  +  ia0R0m-L{,»(^:) 


The  first,  order  vertical  velocity  is  obtained  as  a  solution  of  the  Orr-Sommerfehl  equation 
as  in  Biidges  &  Morris  (1987).  The  metric  m(Y)  is  expanded  in  a  series  of  Ghebyshev 
polynomials, 

m(Y)  =  im„  +  mjT,(F)  +  m2T2(F)  (3.32) 
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wii.Ii  mu  =  m i  =  — |  ami  nt  i  =  L.  Then  the  Chebyshev  pioiluct  ami  differentiation 
formulae  can  be  used  to  write  the  right  hand  side  of  (3.30)  in  a  finite  series  of  Chebyshev 
polynomials.  Similaiily  the  coefficients  ,4(1' )  and  B(Y)  can  be  expanded  in  finite  seres 
of  Chebyshev  polynomials.  It  is  then  straightforward  to  expand  f>2i(K)  in  a  finite 
series.  Substitution  into  (3.30)  and  application  of  the  Chebyshev— r  method  results  in  a 
discrete  system.  After  numerically  eliminating  £  (3.30)  reduces  to  the  Unite  dimensional 
mat  rix  equation 

[D.|(2o0,2/J)J  •  {«2i }  =  {rhs}.  (3.33) 

Assuming  (2tt0,2/f)  is  not  an  eigenvalue  when  (orD,/J)  is  (non-resonance  point)  the  (com¬ 
plex)  system  (3.33)  is  easily  inverted  to  find  the  vector  {O21}  of  Chebyshev  coefficients. 
The  remainder  of  the  second  order  functions  nre  obtained  in  a  similar  fashion.  Then 
the  functions  haj{y)  ave  obtained  by  repeated  use  of  the  Chebyshev  product  and  differ¬ 
entiation  formulae.  The  matrix  (D.j  (rv0,  /?)]  is  the  discretized  Orr-Soimnerfehl  equation 
and  therefore  (when  *«„  is  and  eigenvalue)  has  a  left  eigenvector.  Instead  of  using  the 
adjoint  eigenfunction  of  the  continuous  Orr-Sommerfeld  equation  to  obtain  the  bifurca¬ 
tion  coefficients  we  simply  use  the  left  eigenvector  of  the  matrix  [D.i(a0,/?)].  Let  {>/>*} 
be  the  left  eigenvector  of  D.(;  that  is, 

{<M"-[D.«(cVo,/*)]  =  0,  (3.34) 

then  tli<’  bifurcation  coefficients  h3j  are  easily  obtained  by  taking  the  (discrete)  complex 
inner-prod  net  of  {</>, }  with  [hyj).  Complete  details  of  the  numerical  calculations  will 
be  reported  in  Bridges  [1991b], 

Alternatively  a  complete  analysis  of  transition  process  (through  the  spatial  evo¬ 
lution  equation)  can  be  carried  out  in  the  discrete  setting.  Write  Hie  Navier-Stokes 
equations  as  an  evolution  equation  in  *:  ^4?  =  L(c,  R)  ■  $  -f  N(4';  Jt)  and  for  brevity 
we’ll  sketch  the  2D  case.  Eliminate  the  differential  operators  in  y  by  expanding  4'  is  a 
finite  series  of  Chebyshev  polynomials: 

N 

*(*, y)  =■  5«I’u(*)  +  £  *«(*)r«(y).  (3.35) 

>1=1 


35 


The  evolution  equation  for  $  is  then  reduced  to  an  evolution  equation  on  the  finite 
(although  large)  dimensional  space  R'w, 

~<L>  =  +  N($,/l)  (3.36) 

where  'I'  =  ($oi^i).")5'i  L  is  a  matrix  and  N  is  an  algebraic  nonlinear  operator.  It 
is  then  .straightforward  to  apply  the  usual  bifurcation  theory  for  evolution  equations  in 
finite  dimensions  to  (3.36). 


3.3  Secondary  bifurcations  2D  — »  3D 

An  alternate  route  for  three-dimensionality  to  arise  is  through  a  secondary  bifur¬ 
cation  from  a  2D  state  to  a  3D  state.  In  fact  this  is  a  widely  accepted  theory  for  the 
origin  of  three-dimensionality  (Orszag  &  Patera  [1983],  Bayly,  Orszag  &  Herbert  [1988], 
Herbert  [1988]).  The  "secondary  instability”  theories  of  Orszag  &  Patera  and  Herbert 
are  however  temporal  theories.  In  this  section  we  introduce  a  spatial  theory  for  the 
bifurcation  from  2D  finite  amplitude  states  to  three-dimensional  stales.  The  idea  is  to 
study  the  spatial  Floquet  multipliers  along  a  branch  of  spatially  periodic  2D  states  but 
with  3D  perturbations, 

Let  {v.  v,p)  be  a  2D  spatially  periodic  state  satisfying  (2.1)  and  consider  the  addi¬ 
tion  of  a  31)  perturbation  (u+£,  a  4-  »/,  w,  p+q).  Substitution  into  the?  3D  Navier-Stokes 
equations  (3.1)  and  linearization  about  the  2D  state  results  in  the  following  system  with 
periodic  (in  a:)  coefficients 


°i  +  dJl  +  9A=0 

(lx  dy  dz 


(3.37) 


’«*  +  »&  Uy+Uy 


i]  -f  7  </  -|- 

c/ 


«»+«#  <> 
o  t>4 


(3.38) 

The  pressure  perturbation  q  can  be  eliminated  by  taking  the  divergence  of  the  pertur¬ 
bation  momentum  equations, 


Ary  =  -2  [{Uy  +  Uy)l)r  +  Vyljy  +  +  VX£y}  , 


(3.39) 
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Now,  taking  the  Laplacian  of  the  streatnwisc  and  vertical  velocity  perturbation  equa¬ 
tions  results  in  the  following  coupled  linear  equations  for  (£,»/)  with  periodic  (in  x) 
coelficieuts, 


AA(  -  R-^ -Aq  -  RA  {( Uy  -I-  uy)i)  +  (u  +  U  -  c)£t  +  £ux  + t>£„}  =  0 

A 

AAi)  -  R—Aq  -  RA  [vyj]  -I-  (u  +  U  -  c)qz  +  £vz  +  vi)y}  =  0 


(3.40) 


7 2«>j  +  7  |2^“  ”  ^(u  +  U  ~  «)«’;  j 


'  d'lWj  ,  d2w.  T.  .( 


"1  -  = 
X  Oy  J 


with  At/  given  by  (3.39)  as  a  function  of  (£,»/).  There  is  also  an  additional  decoupled 
eigenvalue  problem  associated  with  the  spatiwise  velocity  perturbation;  in  particular, 
given  (£,jj),  q  is  obtained  from  (3.39)  and  the  spamvise  velocity  perturbation  is  given 
by  I  he  solution  of 

-■sA«.  +  (.  +  t/-c)3j+.^  =  -3j.  (3,11) 

Equations  (3.40)  and  (3.41)  together  form  the  eigenvalue  problem  for  the  secondary 
bifurcations.  But,  since  they  are  decoupled  the  Floquct  eigenvalue  problem  associated 
with  (3.41)  can  be  studied  independently.  Let 

«»(*,} /.-)  =  c7J[«M(*,l/)cosj02  +  u>2(*,It)smj9«), 

then  u>j(x,y)  ( j  =  1,2)  satisfy  the  following  quadratic  in  the  parameter  eigenvalue 
problem, 


0  (3.42) 


(j  ---  1,2)  for  the  spatial  Floquct  exponent  7.  Given  expressions  for  the  periodic  func¬ 
tions  n(x,n)  and  v(x,y)  along  a  branch  of  2D  spatially  periodic  states,  the  eigenvalue 
problem  (3.42)  can  be  discretized  and  solved  using  the  methods  described  in  Section 


The  main  eigenvalue  problem  for  the  spatial  secondary  “instability”  is  the  coupled 
set  (3.40).  Since  the  coefficients  of  (3.40)  are  periodic  in  x  and  independent  of  z  we 


can  take 


teO  - ^  [(£&;{) + (&;’)■ H  (3-43) 
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with  {£} ,  i]j )  ( j  =  1,2)  periodic  in  *.  Substitution  of  (3.43)  into  (3.40)  results  in  the 
nonlinear  (of  degree  4)  in  the  parameter  eigenvalue  problem 

&m)  •  ($;) = £  •  (J;.)  =  o-  (3.44) 

Explicit  expressions  for  the  operators  Lj{e,ft)  are  easily  obtained  by  substituting  (3.43) 
into  (3.39)  and  (3.40).  The  parameter  e  represents  a  parametrizatiou  of  the  branch  of 
2D  stales  and  ft  is  the  spanwisc  wavenumber  of  the  perturbation.  In  particular,  the 
Floquel  exponents  7(e,/J)  will  depend  on  two  parameters  which  will  result  in  potentially 
more  complex  bifurcations  (when  compared  with  the  2D  secondary  bifurcation  problem 
in  Section  2.2  where  ft  =  0).  The  complexity  of  the  bifurcations  can  also  be  increased 
by  changing  parameters  that  alter  the  equilibrium  state  such  as  <•  (the  wavespeed) 
or  using  Falkncr-Skan  flows  rather  than  the  Blasius  flow  for  the  equilibrium  state  or 
adding  a  compliant  wall  along  the  boundary  (Carpenter  &  Gar  rad  [1985]  &  Carpenter 
Si  Mori  is[l 990]). 

The  types  of  secondary  bifurcation  points  to  be  expected  from  the  eigenvalue 
problem  (3.44)  will  be  similar  to  those  described  in  Section  2.2  (  wavelength  doubling 
and  secondary  bifurcation  to  (streamwisc)  quasi-pcriodic  states)  but  with  the  addition 
of  spanwisc  periodicity  of  wavenumber  ft  (which  is  an  independent  parameter). 

To  determine  values  of  the  parameters  (e,  /?)  at  which  secondary  bifurcat  ions  occur 
would  icquire  numerical  solution  of  the  nonlinear  eigenvalue  problem  (3.44).  One  way 
to  show  that  in  fact  secondary  bifurcations  to  spatially  quasi-pcriodic  states  are  to  be 
expected  is  to  introduce  a  second  parameter  whose  variation  brings  the  secondary  bi- 
fm cation  point  down  to  the  origin  forming  a  codimension  2  singularity.  Then  secondary 
bifurcation  to  spatially  quasi-pcriodic  states  can  be  found  in  the  unfolding.  In  Section 
4  it  is  shown  that  for  each  ft  (sufficiently  small)  there  exists  a  codimension  2  point  in 
(c,  R)  space  whose  unfolding  contains  secondary  bifurcations  to  spatially  quasi-pcriodic 
states.  In  fact  all  along  the  upper  branch  of  the  2D  neutral  curve  secondary  bifurcation 
to  quasi-pcriodic  states  (with  periodic  spanwisc  variation)  will  be  prevalent. 

Symmetry  will  play  an  important  role  in  the  secondary  bifurcations.  For  the 
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secondary  bifurcation  to  streamwise  quasi-pcriodic  states  the  normal  form  will  be 
0(2)  x  S’-cqui  variant  with  the  0(2)  action  associated  with  the  spamvise  periodicity 
and  reflection  and  the  S1  action  is  associated  with  the  second  streamwise  wavenum¬ 
ber  (assuming  the  Floquet  multiplier  lies  at  an  irrational  point  on  the  unit  circle).  A 
torus  bifurcation  in  the  continuous  system  is  locally  equivalent  to  a  llopf  bifurcation 
in  a  map  (discrete  time  system).  In  essence  the  above  secondary  bifurcation  is  equiv¬ 
alent.  to  a  Hopf  bifurcation  in  an  0(2)-equivariant  map.  Results  on  llopf  bifurcation 
in  maps  have  been  obtained  by  Chossat  &  Golubitsky  [1988,  p.  1202).  Modulo  some 
(substantial)  technical  details  the  0(2)-cquivaiiant  Hopf  bifurcation  in  maps  resembles 
the  ()(2)-cquivariant.  Hopf  bifurcation  in  continuous  systems.  In  particular  there  will 
be  two  classes  of  secondary  quasi-pcriodic  states.  This  is  easy  to  see  physically:  the 
secondary  slate  can  correspond  to  an  oblique  wave  (where  the  streamwise  wavenumber 
does  not.  resonate  with  the  2D  wavenumber)  or  the  secondary  state  can  be  oriented  in 
the  streamwise  direction  (but  spamvise  periodic  and  again  non-resonant).  The  group- 
theoretic  results  of  Chossat  &  Golubitsky  can  be  used  to  obtain  further  information. 
There  will  be  a  group  orbit,  of  quasi-pcriodic  states;  in  particular  two  oblique  waves 
with  isotropy  subgroup  SO(2)  and  a  continuous  group  orbit  of  “standing”  secondary 
states  (a  toms  of  invariant  tori!)  with  discrete  isotropy  subgroup.  With  the  addition 
of  another  parameter  there  will  also  exist  points  of  tertiary  bifurcation  to  3-tori!  The 
analysis  of  secondary  bifurcation  to  quasi-pcriodic  states  with  symmetry  is  an  interest¬ 
ing  area  for  further  stud}’.  Our  anaiysis  in  Section  4  shows  that  this  bifurcation  will 
play  a  crucial  role  along  the  upper  branch  of  the  2D  neutral  curve  in  shear  flows. 

The  other  class  of  secondary  bifurcations  is  wavelength  doubling  (spatial  Floquet 
multiplier  piissing  through  -1).  Suppose  the  streamwise  wavenumber  of  the  biisic  state 
is  normalized  to  1  (wavelength  2a)  and  that  a  Floquet  multiplier  lies  at  —1.  To  study 
the  bifurcation  of  4^-periodic  states  we  use  the  equations  (3.1)  perturbed  about  the 
2D  state.  On  the  space  of  47r-periodic  functions  however  the  nonlr  ear  problem  is 
ZJ-cquivariant  with  action  p  •  f(  r,y,z)  =  /( x  +  2n,y,z).  The  non-lrivial  spamvise 
variation  however  provides  an  additional  0(2)-equivariance  of  the  nonlinear  equations. 
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Therefore  secondary  bifurcation  to  wavelength  doubling  from  2D  states  to  3D  states 
with  spanwise  periodicity  is  a  Zf?  x  0(2)  equivariant  bifurcation  problem ,  or  in  terms  of 
dynamical  systems  theory,  the  above  bifurcation  corresponds  to  period- doubling  in  the 
presence  of  a  continuous  symmetry.  Period  doubling  with  a  continuous  symmetry  is  a 
difficult  problem.  Acting  on  the  periodic  orbit  with  the  SO(2)  C  0(2)  action  results 
in  a  torus  of  periodic  orbits.  Period-doubling  will  therefore  correspond  to  a  doubling 
of  the  whole  manifold  of  periodic  orbits  (see  Vandcrbauwhedc  [1989,1990]).  In  light  of 
its  importance  for  spatial  bifurcations  in  shear  flows  the  properties  of  period-doubling 
with  0(2)  symmetry  is  an  interesting  area  for  further  research. 

3.4  Secondary  bifurcations  3D  -*  3D 

Secondary  bifurcations  from  3D  states  that  are  periodic  in  the  streamwise  and 
spamvis<  direction  are  of  great  importance  in  shear  flows.  Three  dimensionality  is 
essential  for  true  turbulence  and  the  theory  in  Sections  3.1  and  3.3  presents  two  routes  to 
3D  states  with  spamvise  periodicity.  3D  states  with  spamvisc  periodicity  are  apparently 
an  inevitable  stage  in  the  transition  process.  If  a  (spanwise  periodic)  3D  state  exists 
and  the  st teamwise  flow  is  also  periodic  (it  can  either  be  the  basic  periodic  state  or 
have  wavelength  doubled  (or  n-tuplcd)  any  number  of  times),  then  we  can  apply  two- 
dimensional  spatial  Floquct  theory  to  study  tertiary  bifurcations  in  both  the  streamwise 
and  spanwise  directions. 

Linearization  of  the  set.  of  equations  (3.1)  about  a  3D  state  periodic  in  (*,  z)  results 
in  a  pdc  with  doubly-pcriodic  coefficients  to  which  spatial  Floquct  theory  (Easlham 
[1973,  Chapt.  6])  will  be  applied.  A  basic  question  to  be  addressed  in- the  application  of 
2D  spatial  Floquct  theory  is  the  role  of  spanwise  spatial  bifurcations  in  the  transition 
process.  In  particular,  is  spanwise  periodicity  a  good  assumption  throughout  the  tran¬ 
sition  process  with  bifurcations  essentially  taking  place  in  the  streamwise  direction? 
Or  alternatively,  does  wavelength  n-tupling  and/or  spanwise  quasi-periodicity  play  a 
significant,  lole  in  the  transition  piocess? 
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Suppose  (n,  v,w,2>)  is  a  three-dimensional  state  that  is  periodic  of  period  2n/p  in 
the  span  wise  direction  and  periodic  of  period  27m  (n  €  N)  in  the  streamwise  direction  (it 
doesn’t  inal  ter  how  many  times  the  wavelength  has  doubled  or  n-tuplcd).  Introduction 
of  a  perturbation  about  the  three-dimensional  state  (tt+£,  v+ij,  u>+C>l'+tf)>  substitution 
into  equations  (3.1)  and  linearization  about  the  known  3D  state  results  in  the  following 
linear  pde  with  doubly  periodic  coefficients, 
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(3.46) 

An  equation  for  the  pressure  perturbation  is  obtained  by  taking  the  divergence  of  the 
momentum  equations  resulting  in 


Aq  =  -2 Uyih  -  2 
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Taking  the  Laplacian  of  each  of  the  momentum  equations  in  (3.46)  will  eliminate  the 
pressure  perturbation  via  substitution  of  (3.47)  resulting  in  three  linear  coupled  equa¬ 
tions  with  doubly-periodic  coefficient  for  the  velocity  perturbations  (£,»/,  Q, 


(3.48) 


Two-dimensional  Floquet  theory  (liastham  (1973,  Chapt.  6])  is  easily  applied  to  (3.48). 
Let 

/  €(*.  y>~)  \ 

i)\z,y,:)  =  e71 1+715  f)\x,y,z)  (3.49) 
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where  £, »}  and  (  are  doubly  periodic  of  “period”  2nn  in  x  and  2rr//?  in  z.  Substitution 
of  (3.49)  into  (3.48)  results  in  an  eigenvalue  problem  for  the  pair  of  spatial  Floquet 
exponents  (71,72), 

L(*,-;7.,72)-  =0.  (3.50) 

Equation  (3.50)  is  again  a  nonlinear  in  the  parameter  eigenvalue  problem  in  both  71  and 
72  and  the  degree  is  four.  When  I  he  basic  3D  state  is  known  only  approximately  (as 
say  a  Fouricr-Chebyshcv  series)  the  eigenvalue  problem  (3.50)  will  require  significant 
computational  effort.  However,  the  simpler  question  of  the  role  of  spanwise  bifurcations 
can  be  addressed  by  setting  71  =  0.  Then  (3.50)  is  an  eigenvalue  problem  in  one 
parameter,  the  spanwise  Floquet  exponent  72.  Given  72  as  a  solution  of  (3.50)  with 
71  -=  0  the  spatial  Floquet  multiplier  is  exp[27T72//3]  and  all  the  usual  bifurcations 
(spanwise  wavelength  doubling,  spanwise  quasi-periodicity,  etc.)  are  to  be  expected. 

If  t  he  3D  basic  doubly  periodic  state  is  of  the  standing  variety,  then  Z$  is  in  its 
isotropy  subgroup.  In  other  words  it  is  a  reversible  state  (invariant  under  z  »-»  —2). 
In  this  case  the  spanwise  tertiary  bifurcations  will  be  of  the  type  found  in  reversible 
systems.  The  Floquet  theory  for  reversible  systems  (with  a  reversible  periodic  orbit)  is 
similar  to  Hamiltonian  systems:  periodic  orbits  are  surrounded  by  tori  genetically  and 
n-tupling  bifurcations  (n  >  3)  are  of  codimension  1  (rather  than  of  codimension  2  as 
in  non-teversible  systems).  Consequently  the  bifurcation  structure  of  the  spanwise  ter¬ 
tiary  bifurcations  will  differ  from  the  bifurcations  in  streamwise  direction.  A  spanwise 
Poincare  section  will  have  a  structure  reminiscent  of  a  symplectic  map! 

3.5  Spatial  evolution  of  the  primitive  variables  in  3D 

The  evolution  equation  for  the  primitive  variables  introduced  in  Section  2.3  is  easily 
extended  to  the  three-dimensional  Navicr-Stokcs  equations.  Taking  the  divergence  of 
the  momentum  equations  in  (3.1)  results  in  the  following  Poisson  equation  for  the 
pressure, 


42 


Ap  +  2Uyvx  +  u\  +  v*  +  to ?  +  2(vtuy  +  t utut  +  wyvt)  =  0.  (3.51) 


The  idea  is  l.o  write  the  Poisson  equation  for  p  and  the  spanwise  and  vertical  momentum 
equations  as  evolution  equations  and  use  the  streamwise  momentum  equation  as  a 


constraint.  Let 
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and  the  constraint  induced  by  the  streamwise  momentum  equation  is 
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Assuming  the  functions  (u,  o,  to,/>)  have  periodic  spanwise  variation,  the  evolution  equa¬ 
tion  (3.53)  is  an  0(2)-equivariant  vectorfield  in  a  suitable  function  space. 
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Figure  3.1  Neutral  curves  in  the  (c,  R)  plane  for  the  modified  (real)  3D  Orr-Sommerfeld 
equation  (3.7)  for  /?  >  0. 


4.  Wave  interactions  and  spatially  quasi-periodic  states 


Our  claim  is  that  along  a  branch  of  2D  spatially  periodic  states  secondary  bifurca¬ 
tion  to  spatial  states  that  are  quasi-periodic  in  the  streamwise  direction  (and  periodic 
in  the  spamvise  direction)  are  to  be  expected.  In  this  section  a  theory  is  presented  that 
shows  that  along  the  upper  branch  of  the  2D  neutral  curve  a  secondary  bifurcation  to 
spatially  quasi-periodic  states  is  generic  (in  the  one-parameter  family  of  2D  spatially 
periodic  states).  This  is  shown  by  analyzing  the  codimension  2  singularity  associated 
with  the  inlerscction  of  the  P  =  0  and  p  ^  0  neutral  curves;  in  particular,  the  codimen¬ 
sion  2  singularity  brings  the  secondary  bifurcation  to  quasi-periodic  slates  down  to  the 
origin.  Analysis  of  the  unfolding  of  the  singularity  shows  that  secondary  bifurcation  to 
quasi-periodic  states  is  codimension  1  along  a  branch  of  2D  spatially  periodic  states. 

Suppose  P  is  chosen  arbitrarily  in  the  interval  (0 ,/?mar)  where  pmaT.  ~  .3  (for  the 
Blasius  boundary  layer).  The  neutral  curve  for  >9  =  0  and  ft  ^  0  is  shown  in  Figure 
4.1.  In  the  two-parameter  family  (c, /?)  there  is  a  codimension  2  point,  where  the  two 
neutral  curves  intersect.  At  the  point  (c0,/?0)  the  2D  state  and  the  3D  state  will  have 
pure  imaginary  eigenvalues;  that  is,  at  (c0,i?<,)  the  Orr-Sommerfcld  equation  (3.7)  will 
have  an  eigenvalue  A  =  ia\  (»i  €  R)  when  P  =  0  and  an  eigenvalue  A  =  ia 2  (a2  £  R) 
when  P  £  0.  Note  that  for  each  (i  €  ((),/?,„„*)  there  exists  a  codiinensiou  2  point, 
thcrefoie  such  a  codimension  2  point  can  be  found  at  each  point  on  the  upper  branch 
of  the  2D  neutral  curve.  For  general  p  G  {0,Pmaz)  the  ratio  cri/n2  will  be  irrational 
but  nt  particular  values  of  P  the  ratio  will  be  rational.  We  suppose  henceforth  that 
the  ratio  c»i/a2  is  irrational  and  then  treat  the  codimension  3  points  (c 0,R0,PU)  where 
«i  /cv2  G  Q  as  special  cases. 

As  in  Section  3.1  the  0(2)  symmetry  forces  the  eigenvalue  A  =  m2  (when  p  ^  0) 
to  be  double.  Therefore  the  spatial  centre-manifold  associated  with  the  codimension  2 
point  (c0,l?0)  is  six-dimensional.  What  we  will  show  is  that  the  codimension  2  point 
(c0,  R0)  can  be  treated  as  an  0(2)  equivariant  (spatial)  Hopf-IIopf  mode-interaction 
on  six-dimensions.  The  correspondence  is  useful  because  there  is  an  interesting  normal 
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form  theory  clue  to  Chossat,  Golubitsky  &  Keyfitz  [1986]  (hereafter  CGIv)  that  is 
applicable.  The  theory  of  CGK  is  a  temporal  theory  but  nevertheless  their  existence 
results  will  be  applicable  here  but  the  stability  of  the  bifurcating  qiucsi-periodic  stales 
will  have  to  be  determined  by  other  methods. 


4.1  Bifurcation  of  spatially  quasi-pcriodic  states 

The  evolution  equation  in  the  primitive  variables  given  in  Section  3.5  can  be  re¬ 
cast  as  an  0(2)-equivariant  vcctorficld  in  the  following  way.  Asstiming  periodicity 
(wavenumber/?)  in  the  spamvisc  direction  we  can  write 

CO 

$(*,?/,*)  =  s'l’oOr,!/)  +  ^2  4'2m-i(*,J/)cosm^  +  $2m(a;,j/)sinm/3z  (4.1) 
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The  point,  about  the  (spatial)  Hopf-IIopf  interaction  can  be  made  using  the  linear  part 
of  the  evolution  equation  (in  (3.53)-(3.54)).  Let 
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where  E|$  is  a  G  x  C  matrix  with  unity  at  entry  (4,5)  and  zero  everywhere  else.  With 
the  matrices  Lo  and  !»{{■  the  linear  vectorfield  =  L(c,7l)$  can  b<j  written 
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Let  R.,,,0  =  (  cof sinmflN  ^  ^  ^  ^  |JC  jjie  usuai  rotation  matrix  on 

\-smiiw  cos  m6  ) 

R2.  Then  an  action  for  0(2)  on  the  Fourier  coefficient  space  is  generated  by 

0(2)  =  (n,  K)  (4.7) 

where  7?.  generates  SO(2)  and  is  given  by 

71  =  diag(I0,  J?#  ®  Iu,7l29  ®Io,...,i2m» 
and  1C  generates  Zj?  and  is  given  by 

K.  =  diag(I0,K®Ifl,«®Io,...) 

with  k  --  diag(T, -1).  With  the  action  of  0(2)  given  in  (4.7)  it  is  clear  that  (4.6)  is  an 
0(2)-e(|uivariant  vectorfield  with  0(2)  acting  trivially  on  the  2D  slate  'I’o-  To  study 
the  spatial  eigenvalue  problem  take  =  cXz$>j(y)  then  (4.6)  decouples  into  the 

sequence  of  eigenvalue  problems 
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Suppose  A  6  ff(Lo)  is  a  spatial  Ilopf  bifurcation  point  (Re(A)  =  0  and  Itn(A)  ^  0).  A  t 
(spatial)  Hopf-IIopf  inode-interaction  between  a  2D  state  and  a  3D  state  (equivalently  ' 
a  Hopf-Hoj)f  mode  interaction  with  0(2)  symmetry  on  6-dimensions)  takes  place  if  (for 
the  same  (<:,  7Z)) 
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is  also  a  spatial  Hopf  bifurcation  point  for  some  m.  Without  loss  of  generality  we  can 
take  m  •—  1.  The  sequence  of  eigenvalue  problems  (4.9)  is  equivalent  to  (using  (4.4) 
and  (4.5)) 

«2m  +  (A2  -  »»V2)«2m  =  AT l(U  -  c)Sim  +  Rp2m  (4.11)  j 

lhrn  +  (A2  -  tn2ft2)f> 2m  =  -2\Uyvim 

♦ 

t(>2„,  +  (A2  -  m2ft2)w2m  =  XR(U  -  c) w2m  +  inftRp2m ,  (4.12) 

with  satisfying  the  same  set  of  equations  (i.e.  every  eigenvalue  of  (4.11)-(4.12) 

is  double).  Clearly,  a  mode-interaction  takes  place  (of  six-dimensions)  if  there  exists 
(c0,  T?0)  at  which  (4.11)  has  a  purely  imaginary  eigenvalue  for  both  in  =  0  aud  m  =  1. 
This  is  in  fact  the  case  as  shown  in  Figure  4.1;  that  is,  a  sufficient  condition  is  that  the 
neutral  curves  for  m  =  0  and  m  =  1  have  a  point  of  intersection. 


For  the  bifurcation  at  the  points  (c0,/?0)  where  a  2D  and  3D  state  interact  the 
normal  form  theory  will  be  sketched  with  complete  details  along  with  numerical  evalu¬ 
ation  of  the  coefficients  to  be  reported  in  Bridges  [1991b].  At  the  critical  point  ( c0,Ro ) 
suppose  that  A  =  tori  («i  €  R)  when  ft  =  0  and  X  =  102  («2  6  R)  when  ft  -/-  0.  Then 
the  linear  solution  in  terms  of  the  primitive  variables  at  the  point  (c„,  R0)  is  given  by 
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where  ( ~o  ,“1,-2)  €  C3  are  complex  amplitudes.  Formal  application  of  centre-manifold 
theory  allows  reduction  to  a  vcclorlicld  on  C3.  For  the  normal  form  Proposition  2.3  of 
CGK  is  adapted  to  the  spatial  setting.  At  the  point  (c0,R0)  the  normal  form  is 
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where  ,  J>2,7o,<Ji  and  q>  are  real  functions  of  p,  N  and  A  where  p  =  |j?0|2»  N  = 
|=i|2  +  |=a|2,  A  =  52  and 6  =  |=2p-|=i|2.  Moreover qo (0,0,0)  =  aq  and q\ (0, 0, 0)  =  a 2. 
An  inlet csting  property  of  the  complex  equations  (4.14)  is  that  the  amplitudes  and 
phases  separate:  let  zj  =  rje'+i  j  =  0,1,2  then  (4.14)  decouples  into 


ro=l>o{p)N)A)r0  'j 

l*‘i  =  0>i(/>,JV,A)  +  «p2(p,iV,A))ri  1  (4.15) 

»  J  =  {Vi(p,  Ar,  A)  -  5p2(p,iV,  A))r2  J 


and 


<!>0  =  q0{ptNfA) 

4»i  =  <7i(p,W,  A)  +  Sin 
h  =<li{p,N,A)-  Sq2. 


(4.16) 


The  .set  of  amplitude  equations  in  (4.15)  is  particularity  easy  to  analyze  because  it  is  a 
Z2  Hi  D.,  cquivariant  vectorfield.  The  group  Z2  =  (F0)  and  D*  =  (F|,F2,F)  where 


Fo  •  (»‘o, »’i , »*2)  =  (— r0,ri,r2) 

Fi  ■  (i'0i *‘1  > *’2)  =  (ro,-rj,-r2) 

(4.17) 

F2  •(r0,ri,r2)  =  (r0)r,,  -r2) 

^•(ro.ri.rj)  =  (ro,r2,ri). 

CGK  have  used  group- theoretic  techniques  to  show  the  existence  of  seven  classes  of 
solutions  in  the  normal  form  (4.15).  They  are  listed  in  Table  4.1  along  with  their 
symmet  ry  group  (as  subgroups  of  Z2  ©  D.|).  Types  1,2  and  3  are  the  strictly  periodic 
stales  that,  occur  away  from  the  interaction.  States  5  and  7  require  an  additional 
parameter  (are  codimension  3).  The  interesting  states  are  4  and  G.  They  correspond 
to  two  classes  of  spatially  (strcumwisc)  quasi-pcriodic  states  with  span  wise  periodicity. 
Type  4  involves  interaction  between  an  oblique  travelling  wave  and  a  2D  wave  with 
intependent.  wavenumbers  (but  the.  same  phase  speed)  and  type  G  involves  interaction 
between  a  standing  3D  wave  (actually  travelling  in  the  streamwise  direction  but  periodic 
in  the  spanwi.se  direction)  and  a  2D  wave. 


49 


Further  information  about  the  bifurcating  states  can  be  obtained  from  the  bifur¬ 
cation  equations.  Expansion  of  the  right  hand  side  of  (4.15)  in  a  Taylor  series  and 

truncation  at  order  3  results  in 

Vo  =  coA  +  /*  +  ao  p  +  60Ar 

pi  =  ciA  +  a\p  +  &i  N  (4-18) 

Vi  =  vl 

where  Cii,«o,6o,ci,ai,6j  and  pj  are  real  constants  to  be  determined  and  (A,/t)  arc  the 
unfolding  parameters  of  the  codimension  2  singularity.  CGK  give  a  partial  analysis  of 
the  bifurcation  equations  (4.18).  There  are  numerous  bifurcation  sequences  depending 
on  the  value  of  the  coefficients.  Computation  of  the  coefficients  relevant  to  the  tipper 
branch  of  the  Bhisius  solution  neutral  curve  are  carried  out  in  Bridges  [19911)]. 

Ultimately  the  importance  of  the  spatially  quasi-periodic  states  that  bifurcate 
along  the  upper  branch  of  the  2D  neutral  curve  will  depend  on  whether  they  are  stable 
or  not.  To  determine  stability  of  the  spatial  states  will  require  rcintroductiou  of  time. 
For  the  spatially  periodic  states  (lie  sideband  instability  will  have  to  be  considered 
as  well.  More  generally  the  spatial  states  correspond  to  spatial  invariant  manifolds. 
Therefore  there  will  be  two  steps  in  the  stability  analysis:  stability  with  respect  to 
parametrically  equivalent  manifolds  (i.e.  in  the  spatially  periodic  c;isc,  stability  with 
respect  to  perturbations  of  the  same  wavenumber)  and  secondly  stability  with  respect 
to  other  “nearby”  spatial  invariant  manifolds  (this  is  a  generalization  of  the  sideband 
instability).  A  stability  theory  for  spatially  periodic  state  including  sideband  instability 
is  straightforward  but  a  theory  for  the  stability  of  the  more  complex  spatial  invariant 
manifolds  is  by  no  means  clear  but  is  clearly  of  great  importance  for  determining  the 
stable  spatial  states  in  shear  flows. 


4.2  Spanwise  resonances  and  mode-interactions  on  8-dimensions 

Alt  hough  we  do  not  pursue  it  here  it  is  possible  to  have  interactions  between  two 
3D  stales  resulting  in  an  eight-dimensional  centre-manifold.  Maintaining  the  basic  as¬ 
sumption  of  periodicity  in  the  spanwise  direction  we  look  for  resonances  between  two 
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3D  blocks  (m  >  0)  in  (4.9).  For  example  consider  block  m  =  1  and  m  =  2  which  would 
correspond  to  a  spunwisc  resonance  between  ft  and  2ft.  Figure  4.2  shows  an  example 
of  neutral  curves  in  the  (c,  It)  plane  for  ft  and  2ft  (obtained  by  solving  the  eigenvalue 
problem  (4.11)).  At  the  point  of  intersection  of  the  two  curves  both  ft  and  2ft  result  in  j 
purely  imaginary  eigenvalues  of  (4.11);  that  is,  A  =  i»i  (oq  6  R)  corresponds  to  ft  (at 
fixed  (c, ,,/?„))  and  2ft  results  in  A  =  icv?  (a 2  6  R).  Each  eigenvalue  is  double  due  to 
the  0(2)  symmetry  (and  £*1/02  is  genetically  irrational)  resulting  in  an  8-dimensional 
centre-manifold,  This  interaction  corresponds  to  a  Ilopf-IIopf  mode-interaction  with 
0(2)  symmetry  on  8-dimcnsions  which  has  been  studied  by  Chossat,  Golubitsky  & 
Kcyfilz  (19SG).  The  interaction  is  very  complex  and  produces  quasi-periodic  solutions 
(in  the  present  case  spatially  quasi-petiodic  states)  with  2,  3  and  4  independent  fre¬ 
quencies  (or  wavenumbers  in  the  present  case).  The  8-dimcnsional  mode-interaction 
will  occur  at  higher  Reynold’s  number  than  the  G-dimensional  mode-interaction  and 
therefore  it  would  appear  to  be  of  less  importance.  However,  there  is  an  interesting  dy¬ 
namical  feature  in  the  eight-dimensiona  interaction.  Melbourne,  Chossat  &  Golubitsky 
[19SS]  have  shown  that  hctcradinic  cycles  can  lie  found  in  mode-interactions  on  8-  , 

dimensions.  Attbry,  etal.  [1988)  have  introduced  a  model  for  fully  developed  turbulence  ! 
in  boundary  layers  which  shows  that  heteroclinic  cylclcs  provide  a  good  theoretical  I 
model  for  the  bursting  phenomena.  Therefore  analysis  of  the  8-dimensional  interaction 
in  transitional  boundary  layers  may  provide  a  prelude  to  “dynamical”  behavior  that 
persists  in  fully  developed  turbulent  boundary  layers. 

4.3  Resonant  triads 

In  the  codimension  two  non-resonant  interaction  treated  in  Section  4.1  it  was 
noted  that  at  particular  values  of  ft  the  interaction  is  resonant;  that  is,  there  exists 
distinguished  points  (c 0,n0ift0)  at  which  the  2D  state  has  wavenumber  oq,  the  3D 
state  has  spunwisc  wavenumber  ft„  and  wavenumber  pen/q  where  (p,q)  are  integers. 
The  interesting  resonances  arc  when  p  =  1  and  q  =  2,3  or  4  (strong  resonances). 

It  is  easy  to  show  that  each  of  the  strong  resonances  occur  in  the  Blasius  boundary  ' 
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layer.  Indeed,  the  (/»,</)  =  (1,2)  resonance  corresponds  to  the  Craik  resonant  triad 
although  our  nonlinear  theory  will  differ:  Craik  treats  the  wave  speed  as  complex  and 
shows  that,  the  energy  in  an  unstable  2D  wave  is  transferred  to  the  3D  wave.  In  our 
theory  the  wavespeed  is  taken  as  real  and  the  resonant  interaction  is  treated  as  a  spatial 
(codimcnsion  2  Ilopf)  bifurcation. 

The  resonant,  interactions  an'  demonstrated  as  follows.  Using  the  (c,  11)  neutral 
curve,  pick  a  value  of  o  at  which  there  are  two  ^-intersections  (ll0  and  R\ )  as  shown  in 
Figure  4.3.  Corresponding  to  (c,if«)  is  wavenumber  ou  and  to  (c,/?.|)  is  wavenumber 
a  j .  Now  map  (no,  Ro )  — >  («2, i?i , /l)  with  /?  ^  0  using  the  Squire  transformation.  Then 

n-iRi  =  /?< ion  fv2  =  \Jal  -  ft2. 

Therefore  corresponding  to  Reynolds  number  R\  there  is  a  2D  wave  with  wavenumber 
(Vj  and  a  3D  wave  with  wavenumber  a-j.  Write  their  ratio  as  p  =  a\/ai  then 

_  « i  Ri 

1  (SoRq 

which  is  easily  constructed  from  the  neutral  curve  data  and  is  plotted  as  a  function  of 
c  in  Figure  4.4.  A  resonant  interaction  occurs  whenever  p  =  In  general  p  >  1  but  it 
is  clear  from  Figure  4.4  that  there  exist  values  of  (c,  R,ft)  at  which  p  =  2,3  and  4  but 
in  general  p  will  be  irrational  (corresponding  to  the  states  in  Section  4.1). 

The  normal  form  for  the  resonant,  interactions  p  =  2,3  and  4  will  be  more  difficult 
than  the  normal  form  for  the  non-rcsonant  interaction  (less  symmetry).  The  normal 
form  symmetry  will  be  0(2)  x  S1  on  G-dimcnsions.  Consequently  the  amplitude/phase 
equations  do  not  separate.  Normal  forms  for  resonant  Ilopf-IIopf  interactions  with  0(2) 
symmetry  on  G-dimcnsions  have  not  appeared  in  the  literature.  This  is  an  interesting 
aiea  for  further  research  and  will  lie  of  great  interest  for  understanding  the  How  near 
the  resonant  points  on  the  upper  branch  of  the  2D  neutral  curve. 
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Table  4.1 


Solution  types  in  the  unfolding  of  the 
codimension-2  singularity  (c0,R0) 


Isotropy 
Subgroup  E 

Fix  (E) 

Vectorfleld 
on  Fix  (E) 

Solution  Type 

0 

Z2  x  D4 

0 

_ 

trivial  solution 

1 

{1}  x  D4 

(r0,0,0) 

Po  =  0 

2D  spatially  periodic  state 

2 

Z2x{F2)1} 

(O.ri.O) 

Pi  -  r\p2  =  0 

3D  (oblique)  spatially 
periodic  state 

3 

Z2x{F,1} 

(0,ri,r!) 

Pi  =  0 

3D  (standing)  spatially 
periodic  state 

4 

{1}  x  (F2, 1} 

(ro.H.O) 

po  =  0 

Pi  -  rfp2  =  0 

2D-3D(oblique) 
quasi-periodic  interaction 

5 

z2  x  {1} 

(0,ri,r2) 

Pi  =o 
p2  =  o 

3D(standing)-3D(oblique) 
quasi-periodic  interaction 

6 

{1}  x  {F ,  1} 

(»'o,n,ri) 

•«  3 

►—  O 

II  1! 
o  o 

2D-3D(standing) 
quasi-periodic  interaction 

7 

{1} 

(ro,ri,r2) 

Po  =  0 

pi  =  0 

P2  =  o 

2D-3D(oblique)-3D(standing) 
quasi-periodic  (3-torus)  interaction 

Figure  4.1  Neutral  curve  of  the  Orr-Sommerfeld  equation  for  /?  =  0  and  (J  ^  0 
illustrating  the  codimension  2  intersection  point. 


Figure  4.2  Neutral  curve  for  ft  and  2/5  illustrating  the  interaction  point  for  spanwise 


resonance. 


5.  2D  spatially  quasi-periodic  states  and  the  compliant  wall  problem 


The  introduction  of  a  compliant  wall  to  the  boundary  layer  problem  has  a  two-fold 
effect.  Research  of  Carpenter  &  Garrad  [1985],  Carpenter  &  Morris  [1990]  and  others 
has  demonstrated  that  proper  use  of  a  compliant  wall  (i.e.  optimal  material  properties, 
geometry  and  construction)  leads  to  a  reduction  in  drag  nnd  a  greater  control  over  the 
point  of  transition  in  the  boundary  layer.  On  the  other  hand  the  potential  for  “wall 
dynamics”  leads  to  new  instabilities  and  competing  instabilities  in  the  boundary  layer. 
It  is  this  latter  feature  that  interests  us  here.  A  classification  of  the  instabilities  (and 
the  nonlinear  bifurcations  associated  with  them)  is  of  great  practical  importance  for 
the  design  of  a  compliant  wall  (in  other  words  they  need  to  be  understood  so  they  can 
be  avoided  (or  used  to  advantage!)). 

There  are  a  number  of  singularities  in  the  linear  analysis  of  the  stability  of  the 
compliant  wall  problem  (see  Carpenter  [1990]).  Here  we  consider  a  particular-  singu¬ 
larity  in  the  neutral  curve  for  the  compliant  wall  problem  that  supports  our  theory  of 
secondary  bifurcation  to  spatially  quasi-periodic  states. 

In  Section  4  the  emphasis  was  on  secondary  bifurcation  to  spatially  (streamwise) 
quasi-periodic  state  that  have  spamvise  periodicity.  In  other  words  the  streamwise 
quasi-periodicity  appeared  with  three-dimensionality  and  the  basic  2D  state  was  spa¬ 
tially  periodic.  However,  it  was  shown  in  Section  2.2  that  secondary  bifurcation  from 
spatially  periodic  2D  states  to  spatially  quasi-periodic  2D  states  was  possible.  Identi¬ 
fication  of  the  points  of  secondary  bifurcation  to  quasi-periodic  2D  states  will  require 
numerical  calculation.  An  alternative  is  to  use  the  “codimension-2  strategy”.  The 
secondary  bifurcation  to  2D  quasi-periodic  states  is  of  codimension  1;  that  is,  such 
bifurcations  are  generic  in  the  one-parameter  family  of  2D  spatially  periodic  states. 
In  the  codimension-2  strategy  we  introduce  another  parameter  that  brings  the  sec¬ 
ondary  bifurcation  point  down  to  the  origin.  In  particular  we  will  study  a  singularity 
in  the  neutral  curve  found  by  Carpenter  &  Garrad  [1985,  Figure  11]  by  varying  a  sin¬ 
gle  parameter  (the  elastic  modulus  of  the  wall).  What  we  intend  to  show  is  that  the 
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singularity  found  by  Carpenter  &  Garrad  is  associated  with  a  bifurcation  to  spatially 
quasi-periodic  2D  states.  This  singularity  in  question  may  not  be  of  great  practical 
importance  to  the  compliant  wall  problem  (it  appears  at  a  Reynolds  number  of  about 
3000  (based  on  displacement  thickness)  whereas  the  initial  2D  instability  occurs  at 
R  ~  500)  but  it  is  nevertheless  significant  in  demonstrating  that  bifurcation  to  more 
complex  2D  spatial  structures  (other  than  spatially  periodic)  is  to  be  anticipated. 

Consider  the  linear  coupled  fluid-wall  problem.  The  fluid  is  governed  by  the  linear 
2D  Navicr-Stokes  equations  (in  the  convective  frame)  given  in  Equations  (2.2).  For  a 
rigid  wall  the  wall  boundary  conditions  are  u  =  v  =  0aty  =  0.  For  a  compliant  wall 
the  governing  equation,  when  the  wall  is  modeled  as  a  simple  beam,  is 

Cm'qP  +  +  Cl<E,l  ~  ~Vw  (5.1) 

where  i)(x,  i)  is  the  vertical  displacement  of  the  wall  and  ,  CE  and  C*KE  are  dimen¬ 
sionless  (using  displacement  thickness  variables)  coefficients  and  Pu>(*,0,/)  is  the  fluid 
pressure  at  the  wall.  The  linear  kinematic  conditions  at  the  wall  are  given  by 

v  =  ^  and  «  =  -Uy(0)i].  (5.2) 

With  the  transformation  rwa :  -  c<  the  equations  (5.1)  and  (5.2)  can  be  written  as 


r*  d>v  a.  d2,t  ,  n*  „  -  „  I 

°BdxT  +  c  cMfa2+C,<E  l~~l  w'  J 

The  linear  Navier-Stokcs  equations  (2.2)  together  with  the  boundary  conditions  (5.3) 
(and  appropriate  boundary  conditions  at  infinity)  are  considered  as  an  evolution  equa¬ 
tion  in  x.  Taking  the  dynamical  systems  approach,  let  (u,v,jv/)  =  eAl(ft,i),p,  a/)  then 
the  fluid  equations  reduce  to  the  (modified  (real))  Orr-Sominerfeld  equation, 

+  A2)2  f.  +  XRUyyV  -  \R(U  ~  c)  +  A2)  V  =  0  (5.4) 
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and  the  boundary  conditions  (5.3)  reduce  to 

%  +  =  0 

d3v  dv  A  /A4  \  at  V  =  0.  (5.5) 

df + A2 m + 7  {r' c° + AVCm + CkeR2) 5=0  , 

Equations  (5.4)  and  (5.5)  together  with  appropriate  boundary  conditions  at  infinity 
(v,  v1  -»  0  as  y  — ►  co)  form  a  nonlinear  in  the  parameter  eigenvalue  problem  (of  degree 
5)  for  A  as  a  function  of  c,  R  and  the  wall  parameters  Co,  Cm  and  Cke-  The  variables 
Cb,Cm  and  Cke  are  introduced  to  eliminate  the  streamwise  dependence  of  C'D,  C"M 
and  Cj(B  (see  Carpenter  &  Garrad  [1985,  cqn.  (6.5)]), 

Cj3  =  CbR  3,  C/ce  =  CkeR  and  C m  =  Ca//Z  1 . 


The  boundary  conditions  (5.5)  are  in  a  form  rather  different  from  Carpenter  &  Garrad 
and  Carpenter  &  Morris;  here  we  suppose  c  6  R  and  A  6  C  is  a  spatial  eigenvalue.  In 
Carpenter  &  Morris  the  classic  spatial  stability  approach  is  used;  that  is,  the  frequency 
wSR.ceC  and  a,  the  wavenumber,  is  the  eigenvalue. 

Figure  5.1  shows  a  schematic  of  the  results  of  Carpenter  &  Garrad  [1985]  (taken 
from  their  Figure  11,  p.  498).  By  varying  the  elastic  modulus  (Cm  fixed  and  Cue  and 
Co  varying  dependent  on  E)  the  neutral  curve  varies  dramatically.  In  particular  as  the 
elastic  wall  modulus  is  decreased  there  exists  a  critical  value  of  E  =  E0  at  which  the 
neutral  curve  breaks  into  two  pieces.  Our  interest  is  in  the  critical  point  E  =  E0  when 
the  upper  and  lower  branches  of  the  neutral  curve  first  intersect.  Carpenter  &  Garrad 
treated  the  problem  from  the  temporal  point  of  view  (real  wavenumber  a  with  c  €  C  the 
eigenvalue)  and  plotted  the  neutral  curve  in  a  -  R  space.  In  Figure  5.2  the  singularity 
is  viewed  in  the  (c,  R)  plane  (E0  will  differ  slightly  from  the  temporal  value  in  Figure 
5.1).  In  other  words  there  exist  E  =  E0  at  which  the  upper  and  lower  branch  of  the 
neutral  curve  (in  the  (c,  R)  plane  intersect.  Recall  that  associated  with  each  point  along 
the  neutral  curve  is  an  eigenvalue  A  with  Re(A)  =  0  and  Iin(A)  ^  0.  Therefore  when 
E  =  E0  and  (c,R)  =  (c0,  R0)  the  linear  problem  will  have  eigenvalues  A  =  ±ta i  and 
A  =  dhi«2  with  ai/a2  (generically)  irrational;  that  is,  the  linear  problem  has  spatially 
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quasi-periodic  solutions  and  a  center  subspace  of  4  dimensions.  The  vertical  velocity 
will  have  the  form 


t >(*»!/)  =  2Re  [j45i(y)e'°11  +  (5.6) 

where  A,  B  6  C  are  complex  amplitudes  (note  that  Oj  and  v2  are  in  general  distinct 
functions)  and  the  wall  displacement  will  also  be  spatially  quasi-periodic , 

i](x)  =  -Re  —  At>i(0)cia‘*  +  —Bv2(0)eio’z  (5.7) 

c„  [ai  oc2 

with  related  expressions  for  u(x,y)  and  ;>(z,y). 

The  idea  is  to  apply  centre-manifold  theory  to  reduce  the  spatial  evolution  equation 
to  a  vectorfield  on  R4  to  which  normal  form  theory  is  applied  to  show  the  bifurcation  of 
nonlinear  spatially  quasi-periodic  states.  The  analysis  is  sketched  here  with  complete 
details  to  found  in  Bridges  [1991d].  Let  /  =  (Re(A),Im(/l)>Re(B),Im(5))  6  R4  then 
the  (formal)  centre-manifold  reduction  (as  in  Section  3.1)  can  be  used  to  construct  a 
reduced  vectorfield  for  f; 

A f  =  L(c0,R0).f  +  N(f,Co>R0 )  E  =  E0  (5.7) 


’  0  Orj 
L(c0,n0)=  ~ai  0 


0  a2 
—ot2  0 


The  problem  has  been  reduced  to  a  vectorfield  on  R4  in  which  the  linear  part  has  2 
purely  imaginary  pairs  of  eigenvalues  without  resonance,  a  singularity  that  has  been 
analyzed  by  Takens  and  by  Guckcnheimer  &  Holmes  [1983,  Section  7.5].  If  we  set 
=  |j4|  and  r2  =  |i?|  then  successive  changes  of  variables  reduces  (5.7)  to  a  T2- 
equivariant  normal  form  (to  some  order).  To  third  order  the  normal  form  for  the 
amplitudes  reduces  to 

c!l±  =  n(m  +  r\  +  br*) 

r  (5.8) 

^7  =  r2(ft2  +  cr2  +  dr\)  d  =  ±1 
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(Guckeuheimcr  &  Holmes  equation  (7.5.2))  where  (/ii,/t2)  are  the  two  unfolding  pa¬ 
rameters  (related  to  c  —  c0  and  R  —  R0  in  the  present  case).  Numerical  calculations  of 
the  coefficients  b,c  and  <1  are  carried  out  in  Bridges  [1991d]  for  the  Blasius  boundary 
layer  adjacent  to  a  compliant  wall.  Unfolding  of  the  normal  form  (5.8)  shows  that  the 
two  branches  of  spatially  periodic  states  will  persist.  In  addition  there  are  secondary 
bifurcations  to  quasi-periodic  states  and  if  certain  parametric  conditions  are  met  there 
is  a  tertiary  bifurcation  to  3-tori.  Schematic  bifurcation  diagrams  are  shown  in  Figures 
5.3(a)  and  (b).  In  Figure  5.3(a)  the  secondary  branch  of  quasi-periodic  states  goes  off 
to  infinity  whereas  in  Figure  5.3(b)  the  secondary  branch  connects  the  two  branches 
of  periodic  states  and  includes  a  tertiary  bifurcation  to  a  3-torus.  Note  that  stability 
assignments  are  not  included  in  Figure  5.3.  The  centre-manifold  amplitudes  and  the 
normal  form  in  (5.8)  are  written  in  terms  of  spatial  evolution  so  only  existence  results 
are  obtained.  Determination  of  the  stability  properties  of  the  quasi-periodic  states  is  a 
non-trivial  problem  and  will  require  the  reintroduction  of  time. 

In  obtaining  the  singularity  in  Figure  5.1  Carpenter  &  Garrad  varied  only  one  pa¬ 
rameter.  Note  that  even  in  the  simple  model  of  the  compliant  wall  (5.1)  there  are  three 
independent  parameters  (the  more  sophisticated  models  of  Carpenter  &  Garrad  and 
Carpenter  &  Morris  contain  considerably  more  parameters).  Our  claim  is  that  another 
parameter  can  be  varied  to  bring  the  two  noii'resonant  wavenumbers  in  the  (spatial) 
Hopf  bifurcation  together  as  shown  in  Figure  5.4.  The  configuration  in  Figure  5.4  is 
the  (codimension  3)  1:1  nonsemisimple  Hopf  bifurcation  and  has  been  analyzed  by  van 
Gils,  Krupa  &  Langford  [1990].  This  singularity  is  of  interest  for  two  reasons.  From  a 
practical  point  of  view  the  location  of  high  codimension  singularities  in  the  parameter 
space  of  the  compliant  wall  problem  is  of  interest  in  order  to  “design  around  them”. 
From  a  theoretical  point  of  view  the  1:1  nonsemisimple  (spatial)  Hopf  bifurcation  intro¬ 
duces  new  spatial  bifurcations;  in  particular,  van  Gils,  Krupa  &  Langford  show  that  the 
unfolding  of  the  1:1  nonseinisimple  Hopf  contains  homoclinic  bifurcations  and  period 
doubling  bifurcations  as  well  as  n-tori  (n  =  2  and  3).  Adaptation  to  the  spatial  setting 
will  result  in  interesting  spatial  structures;  in  particular,  the  “spatial  homoclinic”  will 
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correspond  to  a  soliton-like  feature  in  the  shear  flow! 
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Figure  5.2  Neutral  curve  at  the  critical  value  of  the  wall  elastic  modulus  E  — 
the  (c,  R)  plane. 
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Figure  5.3  Schematic  bifurcation  diagrams  for  the  normal  form  in  equation  (5.S) 
showing  how  secondary  bifurcations  to  2-tori  and  3-tori  arise:  (a)  infinite  branch  of  T2 
and  (b)  finite  secondary  branch  of  T2  with  tertiary  bifurcation  to  T3. 
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Figure  5.4  Coalescence  of  the  non-resonant  Hopf-Hopf  interaction  by  the  addition  of 
a  third  parameter  producing  a  1  :  1  non-semisimple  Hopf. 
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